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ABSTRACT 


We present new arcminute-resolution maps of the Cosmic Microwave Background temperature and 
polarization anisotropy from the Atacama Cosmology Telescope, using data taken from 2013-2016 
at 98 and 150 GHz. The maps cover more than 17,000 deg?, the deepest 600 deg? with noise levels 
below 10uK-arcmin. We use the power spectrum derived from almost 6,000 deg? of these maps to 
constrain cosmology. The ACT data enable a measurement of the angular scale of features in both 
the divergence-like polarization and the temperature anisotropy, tracing both the velocity and density 
at last-scattering. From these one can derive the distance to the last-scattering surface and thus infer 
the local expansion rate, Hp. By combining ACT data with large-scale information from WMAP we 
measure Hp = 67.6+1.1km/s/Mpc, at 68% confidence, in excellent agreement with the independently- 
measured Planck satellite estimate (from ACT alone we find Ho — 67.9-£ 1.5 km/s/Mpc). The ACDM 
model provides a good fit to the ACT data, and we find no evidence for deviations: both the spatial 
curvature, and the departure from the standard lensing signal in the spectrum, are zero to within 
1o; the number of relativistic species, the primordial Helium fraction, and the running of the spectral 
index are consistent with ACDM predictions to within 1.5-2c. We compare ACT, WMAP, and Planck 
at the parameter level and find good consistency; we investigate how the constraints on the correlated 
spectral index and baryon density parameters readjust when adding CMB large-scale information that 
ACT does not measure. The DR4 products presented here will be publicly released on the NASA 
Legacy Archive for Microwave Background Data Analysis. 
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1. INTRODUCTION 


The cosmic microwave background (CMB) provides 
the tightest constraints on our model of the universe of 
any current data source. Successive microwave back- 
ground measurements, including WMAP, Planck, the 
South Pole Telescope (SPT), and the Atacama Cosmol- 
ogy Telescope (ACT), have produced increasingly precise 
and detailed maps of the microwave sky in both tempera- 
ture and linear polarization (Bennett et al. 2013; Planck 
Collaboration I 2018; Henning et al. 2018; Louis et al. 
2017). The result has been a remarkably simple stan- 
dard model of cosmology described by only six param- 
eters quantifying the composition of the universe, the 
initial conditions of structure formation, and the reion- 
ization history, most of which are now determined at the 
percent level (Spergel et al. 2003; Planck Collaboration 
VI 2018). 

Despite the simplicity of this model, known as ACDM, 
and the consistency with many other cosmological probes 
(e.g., Cooke et al. 2016; Alam et al. 2017; Abbott et al. 
2018b; Hikage et al. 2019), there are some tensions be- 
tween its inferred parameter values and those estimated 
from low-redshift observations, most notably for the lo- 
cal expansion rate, Ho. While Hy measured from Type 
Ia supernovae distances calibrated using tip-of-the-red- 
giant-branch (TRGB) stars is consistent with the Planck 
estimate (Freedman et al. 2019), the rate measured from 
distances calibrated using Cepheid variable stars (Riess 
et al. 2019), or by using time delays from the HOLiCOW 
strong gravitational lenses (Wong et al. 2019, but see 
recent updates from Birrer et al. 2020 which give differ- 
ent results), is significantly higher than that estimated 
from Planck. Riess et al. (2019) estimates that the devi- 
ation is at the 3-5c level. Determining whether or not a 
new component of the stardard model or new process is 
needed to explain all the current observations is of great 
interest, with new models being explored that could, for 
example, reduce the sound horizon in order to fit Planck 
CMB data with an increased local Hubble expansion rate 
(e.g., Cyr-Racine & Sigurdson 2014; Poulin et al. 2019; 
Knox & Millea 2020). One of ACT’s goals is to make 
a second, experimentally-independent, measurement of 
the inferred Hubble constant with a precision compara- 
ble to Planck to help assess the robustness of the current 
estimate from the high-redshift CMB data. 

The Planck satellite's full-sky maps in many frequency 
bands set the current standard for cosmological precision 
(Planck Collaboration VI 2018). The Planck temper- 
ature map does nearly as well as possible in measuring 
the angular power spectrum of temperature anisotropies: 
cosmic variance dominates measurement uncertainty at 
angular multipoles from £ = 2 to / — 1600. In polariza- 
tion Plancks measurements are signal dominated up to 
£ — 900 in the E-mode polarization power spectrum, and 
up to € = 1150 in the cross-correlation between tempera- 
ture and E-mode polarization (Planck Collaboration XI 
2016). Both SPT and ACT have been extending the 
reach of Planck by making new high resolution, high 
sensitivity temperature and E-mode polarization mea- 
surements (Naess et al. 2014; Louis et al. 2017; Henning 
et al. 2018). From its site in northern Chile, ACT has 
now mapped roughly half the sky. The gain in infor- 
mation beyond Planck is in new temperature measure- 
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ments that constrain damping tail physics at £ > 1500 
as well as secondary signals, and improved E-mode po- 
larization measurements that directly map the velocity 
of the photon-baryon fluid at recombination, providing a 
second view of the physics at last scattering. 

In this paper we present new maps made with ACT 
observations at 98 GHz and 150 GHz during 2013-2016. 
They incorporate data used to make the maps presented 
in Naess et al. (2014) (N14, hereafter) and Louis et al. 
(2017) (L17, hereafter), and additionally include data 
taken during nighttime observations in 2015 and 2016. 
Previous analyses used only 650 deg? of sky; this analy- 
sis produces maps covering 17,000 deg’. We present the 
cosmological parameters derived from these new temper- 
ature and polarization sky maps from ACT, using the 
information contained in their angular power spectra. A 
companion paper, Choi et al. (2020) (C20, hereafter), de- 
tails our methods for estimating and validating the angu- 
lar power spectra and shows the power spectrum results, 
together with such technical considerations as noise char- 
acterization and a wide range of robustness tests both at 
the power spectrum and likelihood level. It also describes 
the procedure to extract foreground-marginalized CMB- 
only spectra starting from the frequency spectra and like- 
lihood. The data products presented in this paper and in 
C20 comprise ACT's fourth data release, DR4, and will 
be available on the NASA Legacy Archive for Microwave 
Background Data Analysis. 

This paper (elsewhere referred to as A20) is organized 
as follows. In 82 and 3 we describe the observations and 
low-level data processing, including ancillary data prod- 
ucts such as beams. 84 presents our current methods for 
making maximum-likelihood maps. The resulting ACT 
DRA maps are unveiled in $5. These maps are then em- 
ployed in C20 to extract CMB power spectra. 96 of this 
paper describes our cosmological CMB-only likelihood 
function and the methods we use to apply it to the power 
spectra from C20. In 87, we show cosmological con- 
straints from ACT data alone and from ACT combined 
with WMAP; these measurements of the CMB sky are 
independent from the Planck measurements. We explore 
extensions to the ACDM model in 88. The concluding 89 
discusses the improvements we can expect from the sub- 
stantially larger ACT data sets collected in 2017-2020. 








2. OVERVIEW OF NIGHTTIME OBSERVATIONS 


ACT is a 6-meter off-axis Gregorian telescope, located 
at an elevation of 5190m on Cerro Toco in the Atacama 
Desert in Chile. In 2013 the ACTPol receiver was first 
deployed, upgrading ACT to a polarization sensitive in- 
strument. Each year from 2013-2015 a new polarization 
array (PA) was deployed. The receiver holds up to three 
PAs, each contained in an optics tube and accompanied 
by a set of three mm-wave silicon optics, a set of corru- 
gated feedhorns, and multiplexing readout for the detec- 
tors. The instrument is described in detail in Thornton 
et al. (2016). 

In this paper we analyse data taken by ACT with the 
ACTPol camera at nighttime, defined to be between 2300 
and 1100 UTC. The observations in 2013 (s13, hereafter), 
initially presented in N14, consist of 648 hours, after data 
cuts, divided between three deep regions D1, D5, and D6 
surveyed with one detector array, PAI, at 150 GHz. In 
2014 (s14, hereafter), observations conducted with two 





Fic. 1.— Mollweide projection of the ACT DR4 footprints in 
Galactic coordinates. The D1, D5, and D6 regions were observed 
in season s13, D56 in seasons s14 and s15, D8 and BN in s15. The 
AA region was observed in s16 and covers roughly 40% of the sky. 
D1 overlaps with BN, D5 and D6 overlap with D56, and all regions 
overlap with AA. Two constant-declination lines are shown in black 
roughly indicating the observable sky from Chile. The grey region 
between Dec — [—60, --22] was not observed, although accessible, 
to maximize overlap with current optical surveys. The grid spacing 
in both Galactic longitude and latitude is 30°. 


150 GHz detector arrays, PA1 and PA2, over a larger 
region called D56 amount to 391 hours after cuts, and 
were initially presented in L17. In 2015 (s15, hereafter) 
a third dichroic detector array, PA3, operating at 98 and 
150 GHz was added. A total of 1271 post-cut hours were 
spent surveying the D56 region, a wider region called BN 
overlapping with the Baryon Oscillation Spectroscopic 
Survey (BOSS) (Albareti et al. 2017), and a deep region 
in the southern equatorial sky called D8. Finally, obser- 
vations in 2016 (s16, hereafter) consisted of 710 hours af- 
ter cuts, spread over a wide region covering roughly 4096 
of the sky, which we call AA (for AdvACT survey). As 
mentioned above, we restrict this analysis to only data 
taken with the ACTPol camera, meaning that for s16 
observations we only consider data from PA2 and PA3. 
Data taken with the dichroic PA4 detector array of the 
AdvACT camera, at 150 and 220 GHz, which replaced 
PA1 in 2016, will be part of a later analysis together with 
data taken from 2017 onwards in the AA region. We also 
neglect the small amount of test data taken with a spin- 
ning half-wave plate installed in front of the PA2 and 
PA3 arrays. 

Table 1 summarizes the ACT s13-s16 observations con- 
sidered in this work, including the number of detectors 
used to make the sky maps, detector array sensitivities, 
and the number of hours of data used in the maps. A 
detailed list of region names and areas can be found in 
95.1. More details regarding scan parameters and survey 
strategy are presented in De Bernardis et al. (2016) and 
reviewed in C20 in the context of null tests. 

Figure 1 shows à Mollweide projection of the ACT 
DR4 region footprints in Galactic coordinates. The new 
ACT observation strategy allows us to independently 
map roughly 50% of the clean sky used for the WMAP 
(Bennett et al. 2013) and Planck (Planck Collaboration 
V 2018) cosmological analyses at higher resolution and 
sensitivity, especially in polarization. In addition, we 
provide CMB data overlapping with on-going and future 
optical surveys, such as BOSS, DESI, DES, HSC, KiDS, 
and Rubin Observatory (Albareti et al. 2017; Dey et al. 
2019; Abbott et al. 2018; Aihara et al. 2019; de Jong, 
Jelte T. A. et al. 2017; Ivezić et al. 2019). 


3. DATA SELECTION AND CHARACTERIZATION 
3.1. Data selection 


The data selection follows the same philosophy used 
for L17, and a detailed description can be found in 
Dünner et al. (2013). Each 10-minute section of the time- 
ordered data from a single detector array (this data unit 
is called TOD) is analyzed in frequency space to charac- 
terize and assess its quality. In the low frequency band 
(~ [0.01—0.1| Hz), the signal is expected to be dominated 
by the atmosphere, which to good approximation is com- 
mon for the whole array. The full correlation matrix is 
computed in the low-frequency band, and detectors that 
do not correlate well with the dominant mode are ex- 
cluded. The response to changes in the focal plane tem- 
perature is the most important contaminant and must be 
handled carefully as it can reach comparable amplitude 
to the atmospheric signal for some TODs. Some detec- 
tors, called *dark detectors," are deliberately not opti- 
cally coupled and used to estimate and deproject thermal 
fluctuations. Additional cuts have been included with re- 
spect to L17. In the high frequency regime (~ [10 — 20] 
Hz), noise properties are computed and detectors are ex- 
cluded based on their white noise level and if significant 
amounts of skewness and kurtosis are detected. 

We mask the brightest sources beforehand in order to 
avoid triggering the glitch-finder or biasing the high fre- 
quency noise analysis. Given the much larger survey area 
compared to L17, we mask 632 sources which are brighter 
than 2mK (roughly 180mJy at 150 GHz). This thresh- 
old is well below the noise level in the TODs. 

Following Diinner et al. (2013), eight quality criteria 
are computed and detectors are excluded if they do not 
meet the requirement for any of the criteria. A detector 
can be flagged for a short period of time or excluded for 
the entire TOD. If more than 40% of the samples for a 
single detector are flagged, the detector is excluded. If 
more than 40% of the detectors are excluded for a given 
TOD, the whole array is discarded. 

The effective number of detectors being mapped varies 
from 60% to 80% depending on the array and season. 
Table 1 shows the median number of detectors used in 
the map-making stage. 


3.2. Relative detector calibration 


The data are calibrated into physical units using the 
equation 


W DA 
ado eek (1) 


where t is a detector index within a single array and fre- 
quency band, t is an index over 10-minute observations, 
qu and di? 2 are the time-ordered data in physical (pico 
watts, pW) and data acquisition units, respectively, Ri 
is the responsivity of the detector measured from a bias 
step operation, f; is a stable optical flatfield factor of or- 
der unity, and gr is an atmospheric correction factor 
of order unity. 

Bias step operations are performed at least every hour, 
and involve a brief modulation of the detector bias volt- 
age while monitoring the detector current. ‘These cali- 
bration operations are analyzed to extract each detec- 
tor’s power-resistance relationship and time constant 
(see next section). The factors Rs; have units of pW 


TABLE 1 
SUMMARY OF S13—S16 OBSERVATIONS, CAMERA, AND INSTRUMENT PERFORMANCE CONSIDERED IN THIS WORK 





Season s13 s14 s15 s16 
Regions D1 D5 D6 | D56 | D56 BN D8 | AA 
PA1 150 GHz 

Median Nget ® 493 475 487 | 456 

Array Sensitivity) [uKys] 16.7 20.2 17.6 

tops ^ [hrs] 114 256 278 

PA2 150 GHz 


Median N get 671 


E Sensitivity [uK vs] 16.5 
tobs [hrs] 912 
PA3 98 GHz 


Median N get 384 


c Sensitivity [uK4/s] 15.9 
tobs [hrs] 704 
PA3 150 GHz 


Median Ndet 
an Sensitivity [uK vs] 
tabs [hrs] 


399 
20.6 
515 


Notes: (a) Median number of detectors considered for map-making. We do not account for the fraction of samples masked within the 
single detector time-stream. In 2015 the addition of the third detector array, PA3, led to an increase of the bath temperature, which 
caused the number of working detectors in PA1 and PA2 to drop. In general, a small loss of detectors is expected from season to season 
due to cryogenic cycles deteriorating the quality of the connections between detectors and readout elements. (b) The array sensitivity 
is estimated on a single 10-minute chunk of time-ordered data (TOD) from the map-maker noise model and inverse-variance averaged 
within the dataset. Variations across seasons for the same array and region are due to fluctuations in atmospheric loading and number of 
mapped detectors. Variations across regions within the same season and array are due to differences in the elevations of the scans (and 
thus atmospheric loading). These numbers can be compared with the ones estimated from planet observation presented in C20. (c) Total 
amount of data mapped, measured in hours of exposure. This number is computed by summing the duration of each data chunk that 
passed TOD-level cuts. It does not consider the fraction of samples within the TOD that were flagged (thus not used), which would lead 


to only a small correction to the number provided. 


per data acquisition unit. The calibration model always 
uses the most recent applicable bias step measurement 
(Grace et al. 2014). 

The factors f; and gm are corrections computed based 


on each detector’s response to the atmospheric common 
mode. Initially, the factors f; and gm are fixed to unity. 


The common mode for each observation is computed as a 
simple mean over detectors, and each detector is then cal- 
ibrated against the common mode, in the low-frequency 


regime, to obtain correction factors po The shape of 


gm at this stage is dominated by a time-invariant con- 


tribution due to the optical efficiency of the detectors. 
This static component is factored out into the flatfield 
factors f;, with mean f; forced to unity. The process of 
computing a common mode and measuring gi?" is then 


iterated once, using the non-trivial flatfield from the first 
iteration. This iteration is performed so that the factors 
atm 


92i" will be (typically) more tightly distributed around a 
value of unity, which is useful for identifying calibration 
problems in individual detectors or observations. In the 
150 GHz data, the f; and iterated gm factors are in- 
cluded in the calibration of data that goes into the map. 
In the 98 GHz data, the flatfield factors f; are included 


in the model, but we continue to fix p = ]. This 


is appropriate because the inferred corrections gr?" ar 


more tightly distributed around unity than is the case 
at 150GHz, except in cases of very good atmospheric 
conditions where the measurements of p themselves 


become quite noisy. 


We note that this process only corrects for relative 
differences in detector calibration. The effects of atmo- 
spheric opacity are included in the absolute calibration 
to planets described in $3.5. 

In L17, detectors for which a responsivity could not be 
extracted from the most recent bias step measurement 
were excluded. We are now using the average respon- 
sivity of the array in such cases, with a correction from 
the atmospheric correlation. ‘This allowed us to recover 
about 2096 of detectors for s14 with respect to L17, with 
a negligible impact on the quality of the relative calibra- 
tion, evaluated using planet observations. 


3.3. Detector time constants 


Our method for characterizing the detector time con- 
stants is summarized in Thornton et al. (2016). Time 
constants are measured using both planet observations 
and bias steps. Time constants computed from planet 
scans, operating using the same optical coupling as the 
CMB scans, are used as baseline values. From these mea- 
surements we find time constants in the range of [1.0—5.4] 
milliseconds (ms) with a median value of 1.8ms. While 
there are only O(100) planet scans throughout the sea- 
son, bias steps, as an alternative, can be acquired more 
often, as frequently as per CMB scan. There is still a dif- 
ference between bias step time constant and planet scan 
time constant, due to the fact that the planet measure- 
ments are a more direct probe of the relevant optical re- 
sponse speed compared with the bias step method. Thus 
a conversion factor is applied to the bias step results to 
account for this. Details are described in Grace et al. 


6 
(2014). 


3.4. Polarization angles 


Following the same procedure used in N14 and L17, we 
compute the polarization angles of the detectors in PAI, 
PA2, and PA3 using a detailed optical model of the tele- 
scope mirrors and lenses in the optical design software 
CODE V! (Koopman et al. 2016; Thornton et al. 2016). 
The angles of individual detectors are lithographically 
defined during fabrication. We propagate the position 
of the detector array, as assembled in the telescope, to 
the sky using a ray trace and compare to observed de- 
tector positions derived from planet observations. This 
determines the physical angle of the detectors within the 
telescope. The optical chain also introduces a polariza- 
tion rotation of up to ~1.7° near the edge of the focal 
plane, furthest from the primary optical axis, which is 
computed in CODE V using a polarization sensitive ray 
trace. This combined with the physical angle of the de- 
tectors in the telescope gives us the final polarization 
angle calibration. We later test for any additional global 
angle offset using the EB power spectrum (see C20 for 
details regarding this test). 


3.5. Calibration with planets 


We obtain the absolute calibration of the instrument 
using measurements of Uranus. This calibration fac- 
tor converts the data from physical units (pW) to CMB 
brightness temperature units? (uK) and is applied after 
the relative detector calibration. The factor is allowed to 
vary between TODs in order to capture temporal fluctu- 
ations in atmospheric transmission but is the same for all 
detectors on an array in the same frequency band. The 
calibration procedure is essentially unchanged with re- 
spect to Hasselfield et al. (2013) and Dünner et al. (2013); 
here we provide a brief summary. 

The absolute calibration factor applied to each TOD 
in frequency band v is modelled as 


T) (2) 


sin Ó 





g,(w,0) — c, exp ( 





Here, w and 0 denote the average precipitable water va- 
por level (PWV) and elevation of the array center, re- 
spectively, during the observation. The values for c, and 
T, are obtained by fitting the model to measurements 
of the absolute calibration taken at different values of 
w/sin@. The individual measurements of absolute cali- 
bration are given by: 


Ty; 
i = ; 3 
J Ape ( | 





where T, denotes the brightness temperature of Uranus 
in uK and A, is the observed amplitude of the planet 
in pW. The index ? runs over Uranus observations. ‘The 
temperature T, is obtained by modeling the planet as a 
disk with uniform brightness temperature that is scaled 


1 https://www.synopsys.com/optical-solutions/ 
codev.html 

? Unless otherwise noted, the units uK refer to CMB differential 
temperature units, where the spectral radiance of a region with 
temperature AT' is S, — (OB,(T)/OT) rog AT, with B,(T) the 
spectral radiance of a blackbody and Temp = 2.725 K the mean 
temperature of the CMB. 


by the ratio of the solid angles of the planet, Qu, and the 
instrumental beam, Qg»: 
Qu 

: 4 
On, (4) 





Tp= Tolv) 


The assumed frequency dependence of Uranus’ bright- 
ness temperature Ty(v) is described in Hasselfield et al. 
(2013). The amplitudes A,; in Eq. (3) are determined 
as the peak amplitude of Uranus in planet maps that 
are made using the maximum likelihood map-making ap- 
proach described in that same paper. 

The parameters c, and 7, in Eq. (2) are fitted per ob- 
serving season, detector array, and frequency band. Sea- 
son s15 and s16 form an exception as they use a shared 
value for T, to improve the model fit. In addition, two 
separate values of c, are fitted for s15: one before and 
one after the telescope was refocused mid season. Finally, 
it should be noted that a correction to the absolute cal- 
ibration is applied at a later stage in the analysis. ‘This 
O(10%) correction is determined from the angular cross- 
correlation between the final ACT maps and the Planck 
HFI maps and is determined per sky region, with overall 
uncertainty of O(1%); see C20 for further details on this 
correction. 


3.6. Pointing corrections 


ACT has a blind pointing accuracy of about 1’ for the 
nighttime observations considered here. This is compara- 
ble to our beam Full-Width at Half Maximum (FWHM) 
of 1.4’ at 150 GHz, so if left uncorrected it would signifi- 
cantly broaden our effective beam. As in L17 we correct 
our array pointing’ by comparing the observed positions 
of bright point sources to their known catalog positions 
for each TOD. However, due to our larger data volume, 
for this analysis we now perform the fit in map-space in- 
stead of in the time domain. The fits are performed by 
building the measured map signal, r = A’ N~'d, and an 
approximation for the inverse of its covariance, C^, for a 
20x20 arcminutes thumbnail around each bright source, 
where d is the TOD data vector, N is the noise matrix 
estimated from the data, and A the pointing matrix (see 
94 for more information about this notation). We then 
minimize 


x! 7 Yn - OP ms y a)" Cri - CL mGs y a) 
(5) 


where i is the index of the source and m(z,y,aj;) is a 
point source model with a pointing offset of (x,y) and a 
peak amplitude of a;. The resulting fit is 10-100 times 
faster than a TOD-level fit, at the cost of slightly lower 
accuracy thus increasing the residual pointing jitter com- 
pared to L17 results (see 83.7). 

Our median single-TOD point source sensitivity is 
30mJy, so only point sources with a flux of 2 150mJy 
contribute usefully to the fit. The likelihood of cover- 
ing a bright point source in a single TOD is proportional 
to the sky area covered by the scan. The typical scan 


3 Detector positions with respect to the array center are well 
measured in the field with planet observations. This method would 
not improve those estimates as point sources are too weak to give 
a reasonable per-detector measurement. 


length^ of 20deg means that about half of the s13-s16 
TODs do not cover a bright source that can be used for 
pointing measurement. We interpolate the missing point- 
ing measurements by fitting a Gaussian process model to 
the available ones. After applying these pointing correc- 
tions we measure a residual per-map pointing jitter and 
include it as a component in our final beam model (see 
93.7 for the details on this procedure). 


3.7. Beams 


We use nighttime observations of Uranus to character- 
ize the instrument’s optical response function, or beam. 
For each observing season and detector array, we re- 
cover both the azimuthally-averaged radial beam pro- 
fles and the corresponding /-dependent window func- 
tions, the latter of which enter directly into the power 
spectrum analysis. While the pipeline by which beams 
are processed is similar to that employed in earlier ACT 
data releases such as N14 and L17, and is more generally 
based on methods described in Hasselfield et al. (2013), 
there are a number of important changes that warrant 
further discussion. 

First, the method by which planet maps are made has 
been significantly improved by changing how large-scale 
correlated noise is treated in the time-domain prior to 
mapping. Instead of modeling and subtracting only a 
single detector common mode from the data (as was done 
for earlier datasets), the new map-making algorithm es- 
timates and removes correlated detector modes within a 
12’ radius of the planet’s position in the same manner 
as the gap-filling technique described in 84.3. This not 
only results in visibly cleaner maps dominated by mostly 
white noise, but it also considerably improves the effec- 
tive signal-to-noise in the outer region (or wing) of the 
beam. ‘The new mapping method does, however, induce 
a non-negligible transfer function in the form of a unique, 
but constant, offset in each map. This offset is estimated 
and removed prior to any subsequent beam analysis. 

There has also been a change in the model used to fit 
the wing of the beam. For a diffraction-limited optical 
system, the wing is expected to decay radially as a power 
law with an exponent of —3, thus motivating a model of 


the form: 
a 


sin? 0 
where the amplitude a and mapping transfer-function off- 
set b are fit parameters. This model does not, however, 
account for the effects of scattering due to the surface 
roughness of the telescope's primary mirror. For an im- 
perfect parabolic reflector with RMS surface deviations 
c, the antenna gain is given by Eq. (9) in Ruze (1966); 
this may be recast in terms of a unit-normalized beam 
B,(@) and corresponding solid angle Q: 


B,(9) a Boe + S(8), (7) 


Q 2TC a 82 82” ; 2 
f cmo P SUNNY — X — — (rc sin 0/ A)* /n 
50 Am ( À ) ° : nn! 2 


4 The ACT s13-s16 dataset is dominated by scans targeting 
regions with an extent of 10-20 deg in the declination direction. 
The s16 AA region uses much longer scans, but they constitute 
only a small part of the data volume. 


B, (0) = +b, (6) 
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where Bo(@) is the unit-normalized beam for an ideal 


reflector, 6? = (47e/X)? is the phase variance, A is the 
effective wavelength of light, c is the surface-error cor- 
relation length, and we have used the relation G(@) = 
(47/Q)B,,(0). While the first term in Eq. (7) is ex- 
pected to decay according to the same power law as a 
diffraction-limited system, the same does not hold true 
for S(0). Using photogrammetry measurements of the 
telescope’s primary mirror, we find that e — 20 wm and 
c — 28 cm, values for which the contribution of scatter- 
ing to the wing cannot be neglected (especially at shorter 
wavelengths). Our new model is thus a modified version 
of Eq. (6) in which the scattering term S(0) has been 
added. 

Finally, there have been a number of modifications to 
how errors are estimated at various points in the beam 
analysis pipeline. When fitting the radial profiles, the 
covariance of the non-linear scaling parameter maz (as 
defined in Hasselfield et al. 2013) with the core linear 
fit coefficients and wing amplitude is now incorporated 
in the overall model uncertainty. This is accomplished 
by using the solution from our original fitting method as 
the starting point for an MCMC sampling of the full pos- 
terior parameter distribution. The resulting chains are 
thinned (to reduce computational load) and transformed 
in the same manner as the radial fit parameters in order 
to recover the full £-space beam transform covariance. 
An independent ‘scattering’ mode is subsequently added 
to this covariance to account for additional uncertainty in 
our new wing model due to possible geometric modeling 
errors’ and temporal fluctuations in the mirror surface 
RMS parameter e in Eq. (8). 

Lastly, as in L17, the covariant errors are further ad- 
justed for uncertainty in the pointing variance correc- 
tions. Unlike previous analyses, however, the contribu- 
tion of position-dependent pointing variance to this un- 
certainty is now estimated by marginalizing over it as a 
parameter (op) in the likelihood: 





L(up, 0p |d, o) = A(0, 0p) e7 E? (Hie)? /2(0? +08), (9) 


where up is the overall pointing variance, A(O, op) = 
IN [2(0? +02)] 7"? is a normalization factor, and {d, ø } 
are pointing variance fits and corresponding errors for N 
bright point sources in the maps. Note that by prop- 
erly including the position dependence in the likelihood, 
we expect this method to yield improved estimates of 
both the total uncertainty and the pointing variance it- 
self. Compared to L17, we find that overall pointing 
variance has increased (e.g., from 9.9 (8.4) to 12.1 (22.5) 
arcsec? in D56 PA1 (PA2)); this may be due to both the 
aforementioned changes as well as lower accuracy TOD- 
level pointing corrections (see §3.6). 

The mean instantaneous /-space window functions for 
each season and array, as well as their corresponding er- 
rors, are shown in Fig. 2. Additional corrections due 
to map-based pointing variance (see Eq. (9)) and RJ- 
to-CMB spectral conversions are applied as in L17. A 
detailed treatment of how the broad passband modifies 
the effective beam depending on source spectral type may 


? We only measure the positions of panel edges in the telescope's 
segmented primary mirror and must therefore model / interpolate 
the deformed surface in order to estimate the true RMS deviation. 
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Fic. 2.— Window functions for the mean instantaneous beam of each array and band in each season. The window functions used for 
interpretation of the survey maps are slightly modified to account for residual pointing variance in the observations contributing to each 
map. The window function errors shown in the bottom panel are strongly correlated between multipoles. 


be found in Madhavacheril et al. (2019) in the context 
of component separation. We release both the mean in- 
stantaneous and map-effective (i.e., corrected) beam ra- 
dial profiles and transforms as part of the current set of 
data products. 


3.8. Polarized Beam Leakage and Sidelobes 


We detect the presence of temperature-to-polarization 
(T-to-P) leakage in both the core of the main beam as 
well as its sidelobes. Although relatively small in mag- 
nitude, the leakage is significant enough that we must 
account for it at various steps in the analysis pipeline. 
For the main beam, we use the same set of Uranus ob- 
servations from §3.7 to measure a polarization response 
in both Q and U Stokes parameters for each detector ar- 
ray and observing season. Since we do not expect Uranus 
to be polarized across our observing bands, we interpret 
these responses as leakages. 

In order to extract the /-space T-to-P leakage func- 
tions, we first convert each set of polarization maps to 
{Qr, Ur} (the third defined basis in Ludwig 1973) via lo- 
cal map-space linear combinations of Q and U. The new 
basis has the convenient property that its azimuthally 
averaged radial profiles have a direct correspondence to 


6 Measurements of Uranus at visible and near-infrared wave- 
lengths have shown evidence of radial polarization patterns in the 
limbs, but a disk-integrated polarization of less than 0.0596 (Schmid 
et al. 2006); while no equivalent data exists to our knowledge, the 
relevant scattering effects in the planetary atmosphere are expected 
to be much weaker at millimeter wavelengths, resulting in net po- 
larization levels substantially lower than this bound. 


the E and B harmonic transforms. In the flat-sky ap- 
proximation, this takes the form of a second-order Hankel 
transform: 


(Eo), B()) — -2 [(Q.6). 0.0) 4006) 049... 0 


Since Q and U are formed by taking linear combinations 
of the individual detector beams within an array, we use 
the same framework to model radial profiles in Q, and 
U, as we do in temperature. This permits, with the 
exception of Eq. (10), the use of identical pipelines for 
processing the core beam response in both temperature 
and polarization. The resulting leakages and their asso- 
ciated covariances are then incorporated in the overall 
power spectrum analysis as described in C20. 

In L17 we described polarized sidelobes of the main 
beam in the PA1 and PA2 arrays. The sidelobes, al- 
though weak in amplitude, were strong enough to cause 
noticeable T-to-P leakage (predominantly with E-mode 
structure). Although L17 did not include any data from 
PA3 detector array, we commented there that polarized 
sidelobes were not detected in PA3. However, with a 
more careful analysis using additional observations of 
Saturn, we confirm the presence of weak polarized side- 
lobes in the 150 GHz detectors of PA3. These sidelobes 
are observed to be approximately one tenth as strong as 
the ones in PA1 and PA2. Maps of the sidelobes are 
shown in Fig. 3, along with the temperature to polariza- 
tion leakage function. 

For the present work, we apply the same sidelobe mit- 
igation technique as was applied in L17, for all 150 GHz 





0.10 
40 
0.05 
_ 20 0.00 
[- 
E 
O 0 —0.05 
5 
> 
—20 —0.10 
—0.15 
—40 
—0.20 


—40 -—20 0 20 40 
X (arcmin) 


: — TE leakage ` 
: TB leakage - 


bi? (1073) 





0 1000 2000 3000 4000 
l 


Fic. 3.— Weak polarized sidelobes in the PA3 150 GHz detectors. 
(Top panel) Map of the polarized sidelobes, obtained by stacking 
20 observations of Saturn. The main lobe of the beam is at the 
origin and is masked in this image. The grayscale is linear in units 
of 1/1000 of the main beam peak, with positive (negative) num- 
bers indicating polarization parallel (perpendicular) to the radial 
direction. The sidelobes consist of four groups of compact features, 
at 30 to 40 arcminutes from the main beam. (Bottom panel) The 
effect of unmitigated sidelobes on CMB spectrum measurements, 
expressed as transfer beam functions on and bo . The side- 
lobes are strongly polarized in a direction perpendicular to the line 
connecting them to the main beam; this results in leakage that is 
primarily from temperature into E-mode polarization. Note that 
these sidelobes are projected out of the data prior to mapping the 
CMB, so the net leakage in the resulting spectra will be strongly 
suppressed relative to what is shown here. 


detectors. This involves modeling the sidelobes as a sum 
of copies of the main beam lobe, and then using that 
model to deproject the sidelobe signal from the time- 
ordered data prior to map-making. 

The PA3 sidelobes have the same general features as 
the ones in PA1 and PA2; they consist of a small group 
of compact lobes, with approximate four-fold symmetry, 
strongly polarized perpendicular to the radius from the 
beam center. For any individual detector, a sidelobe will 
only map to the sky in cases where the main beam is 
within the array field of view. To explain this latter 
behavior more explicitly, imagine a fixed pointing of the 
boresight and consider the function f(x) that gives the 
position on the sky of the camera's peak response, for 
a detector located at position x in the focal plane. The 
signal in detector t originates mainly from the sky near 
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n(z;), but also contains contributions from each compact 
sidelobe 7 at f(x; -- 5;). However, sidelobe j is observed 
for detector 7 only in cases where x; + 8; is within the 
physical bounds of the focal plane. This is consistent 
with an optical effect originating inside the camera rather 
than due to optical elements outside the cryostat. 

The PA3 sidelobes are located at radii of approxi- 
mately 30 to 40 arcminutes, leading to two significant 
differences relative to PA1 and PA2. First, the broader 
stride of the polarized sidelobes means that the angular 
scale of strongest T-to-P leakage is pushed to a lower 
multipole: € = 300 for PA3 as compared with 500 for 
PA] and PA2. Second, because sidelobes only map to 
the sky when the main beam is also in the array field 
of view, the larger radii of the sidelobes leads to fewer 
detectors suffering from each sidelobe. 

We do not detect sidelobes in PA3 98 GHz, even though 
we would have sufficient signal to detect them if they 
had similar amplitude and positions in this band. One 
hypothesis for the sidelobes is that they are a diffraction 
effect induced by grids that make up the various optical 
filters. In such a case, one would expect sidelobes to oc- 
cur at larger radius for lower frequencies. Sidelobes in 
98 GHz would appear starting at a distance of approxi- 
mately 47 arcminutes from the main beam. Since this is 
roughly the size of the focal plane, the sidelobes would 
map to the sky for only a small fraction of detectors, at 
the edges of the focal plane. 


4. MAP-MAKING 


As done for previous ACT data releases and described 
in Diinner et al. (2013), N14, and L17, we adopt the 
maximum-likelihood map-making method to reduce sev- 
eral terabytes of TODs into high signal-to-noise sky 
maps. In summary, this is achieved by solving the map- 
ping equation 


(A* N~1A)m = A* N~'d, (11) 


where A is the generalized pointing matrix projecting 
the data from map to time domain, N is the detector- 
detector noise covariance matrix, which is defined in 
Fourier space and measured from the data, d is the con- 
ditioned and slightly filtered time-ordered data, and m 
is a set of I, Q, and U Stokes parameter maps of the 
sky. In other words, by solving Eq. (11) we find the 
maximum -likelihood solution of Nyix - Nstokes degrees of 
freedom.’ The slight filtering applied to the time-ordered 
data due to ground removal (see §4.4) is not handled in 
a maximum-likelihood framework. The details of the fil- 
tering and the resulting transfer function are discussed 
in the following sections. 

The dataset described in 82 is processed in separate 
bundles in order to produce four (or two for s16) inde- 
pendent split maps for each season, array, frequency, and 
region with uniform noise properties (see $5 for a detailed 
description of the map-based products that are released 
with this work). T'he splitting of observations into inde- 
pendent bundles is the result of an optimization process 
to minimize the noise bias in the power spectrum. TODs 
are grouped into single-day blocks. The distribution of 
such blocks among the 2 or 4 splits is computed by op- 
timizing for (i) even distribution of the data among the 


” The concept of degrees of freedom is further used in 84.5. 
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splits and (ii) uniformity of hit counts across the sky 
for each split. One important constraint imposed on this 
procedure is that observations from different arrays taken 
on the same day belong to the same split (e.g., s14-pal- 
splitO connected to s14-pa2-split0) to confine cross-split 
correlations due to atmospheric noise to a single pair. 

For this analysis, the Enki map-maker,® which was 
developed and tested for N14 and L17 and easily scaled 
to handle large high-resolution maps, is used to produce 
the nominal maps. All previous ACT data releases em- 
ployed the Ninkasi map-maker.? These two software 
codes have been compared and shown to lead to consis- 
tent results at the power spectrum level when used to 
map data in the D56 region. In the following sections we 
describe improvements in the data processing between 
this analysis and L17. 


4.1. Pixelization 


We continue to produce sky maps in equatorial cylin- 
drical projection. We adopt the plate carrée (CAR) pix- 
elization, as opposed to the CEA pixelization which was 
used in previous releases. For CAR, pixels are equally 
spaced in latitude, and have the same number of pixels 
per latitude ring. Although the pixels are no longer equal 
area, CAR avoids the elongation of pixels at extreme de- 
clinations and thus requires fewer total pixels to achieve 
the necessary resolution across the full declination range 
of our maps. At the equator, the resolution continues 
to be 0.5 arcminutes for both 98 and 150 GHz. Exact 
weights for spherical harmonic transforms are provided 
in the libsharp library'?(Reinecke & Seljebotn 2013) 
for the CAR pixelization. 


4.2. Sidelobes and cuts 


As in L17, we estimate the contamination from the 
Moon being observed through the telescope sidelobes 
by mapping the data in Moon-centered coordinates. In 
this case, the data are bundled such that single-season, 
single-array, and single-frequency contamination maps 
are produced. Although the sky signal gradually aver- 
ages away in such maps, we decide to not subtract the 
measured contamination, but rather to fully cut contam- 
inated data. In L17, TODs that were partially contam- 
inated were excluded from the map-making run. In this 
new analysis, individual samples that map onto a region 
of the corresponding moon-centered map with a polariza- 
tion signal > 200 uK are cut. The recovered data were 
shown to not bias the signal power spectra when com- 
pared to previous results. 


4.3. Gap-filling and cuts 


Although the contaminated samples (see also $3.1) are 
not used to find the maximum-likelihood solution map, 
m, in Eq. (11), estimating the noise model, N, requires 
continuity of the detector data. In this analysis, we adopt 
a new method to fill in the missing data samples for noise 


5 GitHub repository: https://github.com/amaurea/enlib 
9 GitHub repositories: https://github.com/sievers/ 


ninkasi c and  https://github.com/sievers/ninkasi 


OCLave 

10 We use the pixell library (https://github.com/ 
simonsobs/pixell) for manipulating CAR maps and per- 
forming curved sky Fourier operations. 


estimation purposes. We compute the detector-detector 
covariance matrix from the non-cut samples and we then 
use it to predict the value of the cut ones. This new 
method improved the estimation of the map-maker noise 
model and reduced the noise power spectra by roughly 
a factor of 1.5 in the 1000 « £ < 2000 region, while still 
leading to unbiased signal spectra. 

In addition, we use this filling method for glitch-like 
(short) cut regions. The original data from longer cut 
regions is kept to estimate the noise model, but still ex- 
cluded when solving for the signal map. 


4.4. Ground subtraction 


To characterize spurious emission from the ground and 
other scan-synchronous signals we adopt a similar analy- 
sis as in L17. Each of the original data splits (4 or 2 splits 
per season, array, frequency and sky region) are fur- 
ther split to separate observations at different elevations. 
For each detector, we then solve for a one-dimensional 
maximum-likelihood map with a one arcminute reso- 
lution along the azimuth coordinate. Any signal that 
changes while performing constant elevation scans, such 
as the CMB, will average away over many observations. 

These single-detector maps, called azimuth pickup 
maps, are used as a template (with scaling equal to unity) 
to subtract the contamination signal before mapping the 
data in the usual sky coordinates. In L17, the azimuth 
pickup maps were used to define data cuts (similarly to 
the treatment of sidelobes) for the D56 region, but the 
same procedure would cause large gradients in the hit 
counts for BN and AA, thus degrading the uniformity 
of the map’s noise properties. Because we now use the 
full information of the azimuth pickup maps, rather than 
simply defining a threshold-based binary mask, we do not 
perform any scan-synchronous cleaning based on polyno- 
mial fitting as adopted in previous ACT data releases. In 
addition, we are able to capture and remove ground con- 
tamination on smaller scales. An aggressive polynominal 
fitting could achieve similar performance, but it would 
heavily bias the maps. 

As mentioned above, it is important to verify that the 
CMB signal is sufficiently suppressed in the pickup maps. 
Any residual CMB signal will be removed from the TODs 
and bias the sky maps. To characterize the transfer func- 
tion induced by this new ground subtraction, as well as 
testing for possible convergence issues and noise bias, we 
inject a simulated I, Q,U map into the TODs and per- 
form an end-to-end mapping run (including estimating 
the pickup maps). For C20, the power spectra of the in- 
put map and the processed map are compared to quantify 
the bias as function of the angular scale, showing a small 
0.376 bias for £ 7 350 (see C20 for more details). 


4.5. Point source treatment 


Small errors in our data model, such as residual point- 
ing jitter, gain errors, intrinsic variability of point sources 
or the implicit assumption that the sky is pixelized, can 
introduce bias in the maximum-likelihood sky maps. For 
ACT, these biases result in O(10~*) of the signal leak- 
ing by a noise correlation length (~ 1deg) away from 
high-contrast regions, typically manifesting as a faint X- 
like pattern around strong point sources. We previously 
handled this by subtracting a point source model from 
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TABLE 2 
SUMMARY OF REGION SIZES AND MAP INVERSE-VARIANCE AVERAGED WHITE-NOISE LEVELS IN 4,IK-ARCMIN 


Season s13 s14 

Array and Frequency PAli5o | PAli50  PA2150 
D1 (131 deg^) 18.6 

D5 (157 deg?) 16.4 

D6 (135 deg?) 12.6 


D56 (834 deg?) 
D8 (248 deg?) 
BN (3157 deg?) 
AA (17044 deg?) 


PA1150 





s15 


s16 
PA2150 PA39g PA3150 PA2 PA39g PA3150 


78.7 118.5 


Notes: The region areas include pixels that (i) have been observed in all seasons and by all arrays of interest and (ii) have an inverse-variance 


weight c 2 0.105%. in each dataset (single-array, single-season, single-frequency) separately. These areas are somewhat bigger (depending 


pix ~ 


on the region) than those reported in C20 and use to compute power spectra. The inverse-variance weighted map noise levels are computed 
within the areas defined above and reported in units of uK-arcmin. Combined map noise levels for each region are provided in the text in 
95. Sub-regions of these areas are selected for the power spectrum analysis in C20. 


the time-ordered data before making the maps, and then 
optionally adding them back again in pixel space. How- 
ever, this procedure only eliminates one of the sources 
of model error: the assumption that the sky is pixelized. 
Due to having more bright point sources, with the addi- 
tion of the BN and AA regions, a new detector array with 
a less symmetric beam (PA3), and a faster but slightly 
less accurate pointing model, we chose a more thorough 
approach this time. 

For point sources with a peak amplitude brighter than 
10mK (roughly 800mJy at 150GHz), we add to the 
maximum-likelihood solution in Eq. (11) an extra degree 
of freedom for each TOD sample where the contribution 
from the source is brighter than 100 uK, which typically 
is true for samples mapping 5 arcminutes around the 
source. By having a degree of freedom for each data 
point, no explicit model is fit to these point sources, and 
thus the uncertainties associated to the beam model, 
the amplitude of the source as function of time, and 
the pointing model do not introduce any error into our 
treatment. This procedure is applied to 64 (86) sources 
at 98 (150) GHz, corresponding to roughly 0.01% of the 
sky. In theory, the extra degrees of freedom come at 
the cost of worse handling of correlated noise in such re- 
gions. However, since the noise is approximately white 
on their length scale there is in practice no impact on our 
map noise properties. See Næss (2019) for a more thor- 
ough discussion of model errors and biases in maximum- 
likelihood map-making, and the section on source sub- 
sampling in particular for the method used here. 

For dimmer point sources, but still detected with a sig- 
nal to noise ratio SNR>5 in each per-season, per-array 
map, we perform the point-source model subtraction at 
the TOD level, as mentioned above and described in L17 
and Diinner et al. (2013). This method only removes the 
season-averaged point sources signal. Although this ap- 
proximation is sufficient to guarantee an unbiased map- 
making procedure, further map-space refinements are 
adopted as described in §5.1, Neess et al. (2020), and 
Madhavacheril et al. (2019). 


4.6. Noise iterations 


Because of the wide sky area covered over four differ- 
ent seasons, especially s16, and the increased large-scale 
sensitivity due to integration time, we require stricter 
control of the biases induced by the map-making proce- 


dure. The estimation of the noise covariance matrix, N, 
using the data d introduces a small bias at intermedi- 
ate and large scales. This bias is iteratively removed, as 
described in Diinner et al. (2013), by multi-pass map- 
making where corrections to the map estimate m; from 
pass i are computed by running noise estimation and 
mapping on modified time-ordered data d — Am;. In this 
new analysis as a precaution measure, we perform three 
iterations to correct for such a bias as opposed to two as 
performed in L17. 


5. DRA MAP PRODUCTS 


ACT DR4 maps provide the first signal-dominated 
arcminute-resolution E-mode fluctuation maps over 
thousands of deg? and additional small-scale tempera- 
ture anisotropies over 40% of the sky. The products are 
described in this section and will be publicly released on 
the NASA Legacy Archive Microwave Background Data 
Analysis (LAMBDA) website.'! iPython notebooks de- 
veloped to produce tables and figures in this section can 
be found on the ACT GitHub.'? 


5.1. Products and definitions 


The map-making procedure, described in $4, requires 
roughly 8.3 core-years on a modern high-performance 
computing cluster to reduce 11.3 TB of time-ordered 
data into roughly 300 GB of science-quality sky maps. 
The resulting maps have observed-pixel counts ranging 
roughly from 1.9 million (for D1) to 278 million"? (for 
AA). 

The released maps can be grouped into 20 different 
sets, one for each region, season, array, and frequency 
combination (see Table 2 for details). Each set contains 
4 (2 for AA) separately-processed independent splits and 
one map-based co-add for a total of 94 data units. For 
each unit we release 4 FITS files that fully describe the 
data and its properties: 


e Source-free maps: I (or T), Q, and U Stokes com- 
ponents of the sky in units of microkelvin. The 


ll nttps://lambda.gsfc.nasa.gov/product/act/ 

12 httpsi/7github.com/ACTCollaboration/act drà4 em 

15 Each map has an extra 30-4096 overhead of unobserved (zero) 
pixels due to the scan-strategy and due to RA-Dec region center 
offsets between different detector arrays. 
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FIG. 4.— ACT DR4 survey depth and sky coverage maps for 98 GHz (top panel) and 150 GHz (bottom panel) in Equatorial coordinates. 
Lines at constant map-noise are plotted ranging from 7.5 uK-arcmin(red) to 101 uK-arcmin(blue). The Planck 353 GHz intensity map is 
shown in the background. The x-axis (y-axis) shows the RA (Dec) coordinates in degree. 


polarization components are defined by the IAU 
convention!^ (Hamaker & Bregman 1996).  AI- 
though the signal from point sources is removed 
at the TOD level to avoid biases in the maximum- 
likelihood maps (see $4.5), a second more accurate 
map-space point source subtraction is adopted to 
reduce artifacts that could be induced by map-level 
Fourier filtering often used in cosmological analy- 
ses (e.g., C20; Darwish et al. 2020; Madhavacheril 
et al. 2019; Han et al. 2020). Optimally co-added 
single-frequency maps (Naess et al. 2020), made 
from nighttime ACT data from 2013 to 2017, are 
used to measure the positions of point sources. For 
each source detected with a SNR>5, a symmet- 
ric beam-shaped profile at the measured location 
is used as map-space template. The amplitude of 
the template is fit for in each split map both in tem- 
perature and polarization. This procedure better 
captures small variations in the sources amplitudes 
in each split (see Næss et al. 2020 for more details). 


Source maps: I (or T), Q, and U Stokes compo- 
nents of the point-source signal subtracted follow- 
ing the procedure described above in units of mi- 
crokelvin. The simple sum of the source-free and 
source maps represent the observed sky. 


Inverse-variance map: temperature-temperature 
(or IT) component of the 3 x 3 per-pixel inverse co- 
variance matrix in units of UK ?. The 3 x 3 matrix 
is computed by projecting the diagonal component 
of the map-maker noise model onto the sky (thus 
dominated by the high-frequency part of the noise 


14 This is explicitly specified in the FITS header by the keyword 
POLCCONV-IAU. The ACT DR3 maps were released following the 
COSMO convention (Górski et al. 2005a). 


spectrum). Although the inverse-variance map is 
roughly proportional to a hit count map, it fur- 
ther takes into account the variability in detector 
performances and is thus better suited for map- 
based weighting and map co-addition procedures. 
Because of our scan-strategy and large number of 
detectors, the off-diagonal elements (JQ, IU, and 
QU ) are negligible. In addition, the QQ and UU el- 
ements are each well approximated as having half 
the inverse variance of JJ. For these reasons we 
only release the JJ component of the full 3x38 co- 
variance matrix. This product is used in C20 to 
exclude noisy regions and to weight the spectra. 


Cross-linking maps: T-, Q-, U-like components 
used to describe the cross-linking properties of a 
dataset. These components are computed by pro- 
jecting the time-domain vector, d = const, onto the 
sky using a projection matrix constructed by artifi- 
cially setting all detector angles to zero (i.e parallel 
to the azimuthal motion of the telescope). The T- 
like component is the usual inverse-variance map 
(see above). For a small-throw perfectly azimuthal 
scan over a setting (rising) region of the sky around 
Dec = 0, the resulting footprint in equatorial coor- 
dinates is a nearly-straight scan tilted by roughly 
+45 (—45) degrees with respect to the equator. In 
such configuration the observed pixels in the Q, U- 
like maps will have Q ~ 0 and U ~ +I for the 
setting scan and Q ~ 0 and U - —I for the rising 
scan. In other words, the Q component of the cross- 
linking map mainly encodes the local angle between 
the scanning direction and north, while the U com- 
ponent mainly encodes the balance between rising 
and setting scans. These maps are used in C20 
to select regions with uniform cross-linking, which 


translates into uniform noise properties. 


Although the time-ordered data are calibrated before 
entering the map-making step, an additional calibration 
factor is estimated for each map via a spectrum-level 
comparison with the Planck temperature maps (see 83.5), 
as described in C20. The released maps have been recal- 
ibrated, but no correction for polarization efficiency has 
been applied. Similarly, the power spectra of these maps 
entering the likelihood analysis have been corrected for 
the calibration factors and the covariance matrix includes 
calibration uncertainties. Polarization efficiency is then 
applied correcting the model within the likelihood code 
(see C20). 


5.2. Sky maps 


Figure 4 shows the s13-s16 survey area and the tem- 
perature map noise levels for both the 98 and 150 GHz 
channels. Although observations at 98 GHz only started 
in 2015, the spatial distribution of the noise is similar be- 
tween the two frequencies, with the exception of the D1, 
D5, and D6 regions observed only in s13.!? The Planck 
353 GHz intensity map is also shown, highlighting that 
with the AA survey we have observed a considerable frac- 
tion of the Galactic plane. The Galactic center is part of 
the survey but sparsely sampled in time. 

More quantitatively, Fig. 5 shows the cumulative distri- 
bution function of the map white noise. ‘The information 
about the spatial distribution of the different regions is 
lost in this representation; however it is a useful tool for 
forecasting studies. When combining 98 and 150 GHz 
data, roughly 600 deg? of the sky have been observed at 
a noise level < 10 uK-arcmin. Approximately 4000 deg? 
of the s13-s16 150 GHz survey have a noise level that is 
lower than the full-mission whole-sky average 143 GHz 
Planck channel. The corresponding number for ACT 
98 GHz compared to Planck 100 GHz is 9000 deg?. It is 
worth noting that (i) some of the deepest regions in the 
Planck survey are observed by ACT, (ii) at large angular 
separations the ACT temperature noise strongly departs 
from a white spectrum due to atmospheric correlations; 
much less so for the polarization channels, and (iii) the 
comparison does not account for differences in the beam 
size, which significantly increase the Planck uncertainties 
at small scales compared to ACT. A comparison of the 
ACT and Planck uncertainties at the power spectrum 
level is presented in C20. 

The list of effective region sizes, computed by exclud- 
ing noisy parts at the edge of each region, and inverse- 
variance averaged white noise levels are presented in Ta- 
ble 2. D56 is the best cross-linked region with a map- 
noise level!® of 11.5 (18.4) uK-arcmin at 150 (98) GHz 
over 834deg?. Figures 6 and 7 show a slightly-filtered 
cutout of the D56 region observed at 98 and 150 GHz, 
respectively. The top three panels show the measure- 
ment of the three Stokes parameters T', Q, and U (from 
top to bottom) highlighting scales between 1 degree to 
a few arcminutes. The temperature fields look identical, 
as expected from a high signal-to-noise measurement of 
a low-foreground part of the sky. The Q and U mea- 
surements are also clearly signal dominated and highly 


15 D5 was also observed in s16 and mapped together with the 
AA wide scans for simplicity. 
16 We ignore the contribution of the D5 and D6 regions. 
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Fic. 5.— Map noise cumulative distribution function for 


98 and 150 GHz and their combination. Two horizontal lines show 
the Planck full-sky averaged white-noise level for the 100 GHz and 
143 GHz channels. 


correlated between frequencies with with plus-like and 
cross-like patterns for Q and U that are characteristic 
of E-mode polarization. The two bottom panels show a 
combination of the Q and U fields into E and B fields. 
The E field is clearly signal dominated, at a signal-to- 
noise much higher than previously reported in L17 and 
with the addition of the 98 GHz channel. The B field 
map is clearly noise dominated, thus providing a proxy 
to understand the noise properties of the dataset. The 
D5 (in the top-right corner) and the D6 (in the bottom- 
left corner) regions are now barely visible in the co-added 
map compared to L17. ‘This shows that the s15 data 
have greatly improved the overall depth of the D56 re- 
gion. In addition, a localized 45°-tilted quadrupole is 
clearly visible at RA = 1.56 and Dec = —6.39, with a 
corresponding 0°-tilted quadrupole in the E map, show- 
ing the polarized emission of the bright active galactic 
nucleus (AGN) PKS0003-066. 

The BN region has an extent of 3157 deg? with a noise 
level of 29.2 (33.9) uKK-arcminat 150 (98) GHz, respec- 
tively. The degree of cross-linking of this region is bi- 
modal: roughly 1800 deg?, between right ascensions of 
170 and 235 (or —125) deg, have been observed both dur- 
ing rising and setting, i.e. with the sky at two different 
orientations with respect to the constant-elevation scans 
(see Fig. 4). This subset of BN is used in C20, Mad- 
havacheril et al. (2019), Namikawa et al. (2020) and Dar- 
wish et al. (2020). The remaining 40% of the region was 
observed only during rising scans, and is currently used 
for point-source and cluster searches. Figures 8 and 9 
show a slightly-filtered cutout of the 7’, E, and B fields 
measured in the BN region at 98 and 150 GHz, respec- 
tively. Both frequencies show a signal-dominated mea- 
surement of the E field. 

The AA region covers roughly 17,000deg? of the 
sky reaching a noise level of 62.1 (78.7) uK-arcmin at 
150 (98) GHz. The degree of cross-linking varies substan- 
tially across the observed area. Constant-elevation scans 
from the ACT site project to arcs on the equatorial sky, 
and near declination —35°, these arcs tend to run nearly 
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Fic. 6.— Rectangular cutout of the s13-s16 98 GHz maps in the D56 region. The selected area covers 630 deg? of the sky. To visually 
highlight the wide range of measured angular scales, modes with |/| $150 and |€;|<5 have been filtered out. The top panel shows 


temperature fluctuations in a range of «250 UK. The remaining four black-and-white panels show (top to bottom) Q, U, E-mode, and 
B-mode polarization measurements in a range of +30 uK. The x-axis (y-axis) shows the RA (Dec) coordinates in degree. 
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Fic. 7.— Rectangular cutout of the s13-s16 150 GHz maps in the D56 region. The selected area covers 630deg? of the sky. To 
visually highlight the wide range of measured angular scales, modes with |/| 2 150 and |£;| $ 5 have been filtered out. The top panel shows 


temperature fluctuations in a range of +250 uK. The remaining four black-and-white panels show (top to bottom) Q, U, E-mode, and 
B-mode polarization measurements in a range of 430 uK. The x-axis (y-axis) shows the RA (Dec) coordinates in degree. 
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Fic. 8.— Rectangular cutout of the s13-s16 98 GHz maps in the BN region. The selected area covers 1230 deg? of the sky. To visually 
highlight the wide range of measured angular scales, modes with |/| $150 and |€;|<55 have been filtered out. The top panel shows 


temperature fluctuations in a range of +250 uK. The remaining two black-and-white panels show (top to bottom) E-mode and B-mode 
polarization measurements in a range of +30 uK. The x-axis (y-axis) shows the RA (Dec) coordinates in degree. 


parallel to lines of constant RA, making it difficult to mation in the s13-s16 ACT temperature and polarization 
achieve strong cross-linking. Some variety in trajectory measurements. ‘The left-hand panel shows the emission 
is achieved by observing at different elevations, but the of the five highest signal-to-noise point sources detected 
net result is that mid-latitudes have somewhat differ- at 150 GHz ordered (top to bottom) by decreasing unpo- 
ent noise properties than the higher and lower latitudes larized flux. The size difference, especially visible in the 
where the cross-linking angle is closer to 90 degrees. In temperature thumbnails, is due to the different resolution 
general, it is possible to define stripes at constant declina- at 98 and 150 GHz (see $3.7 for details). The polariza- 
tion in the AA maps that have uniform noise properties tion signal, in Q and U, is also clearly correlated between 
(modulo differences due to exposure time). This feature frequencies. The polarized flux from bright AGNs was 
has been used in the power spectrum analysis presented previously measured using the ACT DR3 150 GHz data 
in C20. and showed to be on average roughly 3% of the unpo- 


The D8 region has comparable map noise levels to the larized flux (Datta et al. 2019). The right-hand panel 
D56 region but over a smaller region of 248 deg? of the shows the five highest signal-to-noise SZ-detected clus- 
sky. Given its location on the sky, the degree of cross- ters which, as expected from the tSZ effect, produce a 
linking is very low. temperature decrement both at 98 and 150GHz. A de- 

Figure 10 illustrates some of the high-resolution infor- tailed analysis of SZ-selected clusters will be presented 
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Fic. 9.— Rectangular cutout of the s13-s16 150 GHz maps in the BN region. The selected area covers 1230 deg? of the sky. To visually 
highlight the wide range of measured angular scales, modes with |@| 5150 and |€;|<55 have been filtered out. The top panel shows 


temperature fluctuations in a range of +250 uK. The remaining two black-and-white panels show (top to bottom) E-mode and B-mode 
polarization measurements in a range of +30 uK. The x-axis (y-axis) shows the RA (Dec) coordinates in degree. 


in Hilton et al. (2020). filtered using a radially-symmetric beam as signal tem- 

In summary, the combination of the ACT DR4 map plate!’ and a 2D correlated noise model measured from 
products provides the first signal-dominated arcminute- the same map. We then look for positive temperature 
resolution temperature and E-mode fluctuation maps fluctuations with a detection threshold of SNR > 2 (2.8 
over 4000 deg? of the sky. They also add substantial for AA). At 150 GHz, each region of the full s13-s16 
small-scale temperature information over 40% of the survey is observed by at least 2 independent datasets 
sky, complementing existing Planck measurements. (maps). We use this feature to clean the single-map cat- 


alogs from random noise fluctuations by cross-matching 
all the 16 catalogs and requiring a source to be observed 
5.3. Point source masks at least twice. A final catalog of source positions, SNR, 
and fluxes is then formed by inverse-variance weighting 


Alth h i | f h hvsical | if ae 
though interesting trom the astrophysical point o Bend dudas loc 


view, the presence of bright point sources in the maps 
ee ener of bias for many cosmological analyses. For 17 Each map is paired with the corresponding jitter-corrected 
this reason we release a set of point source masks as part beam (see 83.7). 

of the DR4 data products. Each of the 16 per-region, 18 The weighting is done by assuming independent errors. This 
single-season, single-array maps at 150 GHz is matched is a good approximation at arcminute scales where the noise is 
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Fic. 10.— Sample of compact arcminute-scale objects detected in the ACT DR4 maps. All thumbnails are extracted by reprojecting a 
subset of the s13-s16 single-frequency map onto the plane tangent to the source coordinates to remove declination-dependent distortions. 
Each thumbnail has an extent of 15 x 15 arcminutes and a resolution of 0.larcmin. (Left panel) Top-five highest signal-to-noise point 
sources detected in temperature at 150 GHz, ordered (top to bottom) by decreasing peak temperature. The colors (white to black) cover 
a range of [0,5] x 10^,K in temperature and +1.5 x 10°yK in polarization. No correction for the T-to-P leakage from the main-beam 
was applied. (Right panel) Top-five highest signal-to-noise SZ-selected galaxy clusters detected in a combination of the 98 and 150 GHz 
temperature maps, ordered (top to bottom) by decreasing signal-to-noise ratio. The colors (white to black) cover a range of [—500, 250] uK. 
An offset computed for each thumbnail (excluding the central source) is removed. 


The catalog is effectively complete for sources detected 
at a combined SNR 5, but because of the varying map 
depth, the corresponding flux limit is a function of sky 
position. To take full advantage of the deep regions and 
to keep the number of masks at a minimum, we produce 
two sets of masks: (i) one set with all sources with flux 
> 100mJy over the whole AA footprint including BN, 
and (ii) another set with all sources with a flux > 15 mJy 
only covering the D1, D56, and D8 regions. Both sets are 
produced for 5, 8, 10, 35 arcmin hole radii. The first set 
contains 1215 holes over roughly 17,000 deg?, the second 
one 1012 holes over roughly 1,200 deg?. 

A visual inspection of the DRA maps shows the pres- 
ence of a few extended astrophysical objects!?. In some 
cases the center of the objects is detected by the point 
source finder, but the extended structures may not be 
fully masked. In addition, some objects are too faint to 
be robustly detected. We therefore additionally include 
an extended-source mask. Using the HyperLeda catalog 
(Makarov et al. 2014), we select 414 objects with a mean 
surface brightness «25 mag in the optical B-band and 
with a linear extent d4'.^? Extended objects identified 
by this procedure include members of the Virgo cluster. 
We then mask a radius d/2 + 2 arcmin at the source lo- 
cation. Both versions of the masks are used as input in 
the C20 power spectrum analysis. 


greater than the residual CMB signal. 

19 A description of such objects can be found in Ness et al. 
(2020). 

20 We pick as linear extent the maximum between the RA and 
Dec angular extent provided in the catalog. 


6. LIKELIHOOD AND PARAMETER ESTIMATION 
METHODS 


6.1. ACT DRZ power spectra and likelihood 


C20 describes the estimation of the power spectra of 
the DR4 maps and their covariance matrix. Table 3 
summarizes the regions, maps and frequencies that are 
used, and the number of individual cross-spectra that en- 
ter each region’s power spectrum. As described in C20, 
the spectra for regions which have sources masked to a 
15 mJy level are co-added into a ‘deep’ product, with the 
‘wide’ product composed of spectra from regions masked 
to a 100 mJy level. The foreground-marginalized CMB 
bandpowers (in TT, TE, EE) are then estimated from 
the multi-frequency spectra, including data from ACT- 
MBAC DR2 (Das et al. 2014), for these two subsets of 
the data. The resultant DR4 coadded spectra are shown 
in Fig. 24 of C20 and in 87 of this paper. 

The likelihood of the derived CMB-only ACT band- 
powers can be well described by a simple Gaussian dis- 
tribution, summing the contribution from the deep and 
wide bandpowers as 


—2InLacr = —2InLact,d — A2lnCACT,w ; (12) 
with 
—2InZ4cr,a/« — (Ct, — gy EP - O7), 
(13) 
where Oe is the binned CMB theory, C /* are the deep 
and wide CMB bandpowers and %iq/,, are their covari- 


ance matrices. Although the deep and wide spectra 
are both estimates of the same underlying cleaned CMB 
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TABLE 3 
REGIONS, AREAS, MAPS AND FREQUENCIES THAT ARE USED IN THE DR4 POWER SPECTRA (DESCRIBED IN C20), AND THE NUMBER OF 
INDIVIDUAL CROSS-SPECTRA THAT ENTER EACH REGION’S POWER SPECTRUM. 


Region Area Maps? 
(deg?) 

Deep* D1 23 s13 PAI 
D5 20 s13 PA1 
D6 20 s13 PA1 
D56 340 s14 PA12, s15 PA123 

Wide BN 1400 s15 PA123 
AA-w0 570 s16 PA23 
AA-wl 1170 s16 PA23 
AA-w3 880 s16 PA23 
AA-w4 1070 s16 PA23 
AA-w5 210 s16 PA23 


Timaps 
98 GHz 





98 x 98 98 x 150 150 x 150 
150GHz | Nrr® Nre Nrr Nre Nrr Nre 
0 0 0 0 1 1 
0 0 0 0 1 1 
0 0 0 0 1 1 
1 1 5 10 15 25 
1 1 3 6 6 9 
1 1 2 4 3 4 
1 1 2 4 3 4 
1 1 2 4 3 4 
1 1 2 4 3 4 
1 1 2 4 3 4 


a The maps are labeled by season of observation (s13, s14, s15 or s16) and by detector array (PA1, PA2, PA3). When multiple arrays 
are used, PA123 indicates PA1+PA2+PA3, for example. Only PA3 has detectors at 98 GHz. The number of maps is nmaps. 


^ 'The number of TT cross-spectra, NT, is equal to Ngg. 


© The ‘deep’ regions have a point source mask with a 15 mJy threshold; the ‘wide’ regions have a 100 mJy threshold mask. 


power spectrum, their bandpower window functions that 
bin the theory are different (see C20), and so are kept 
separate in the likelihood. 

The covariance of the CMB bandpowers already ac- 
counts for foreground uncertainty, and includes calibra- 
tion and beam uncertainty. Only one nuisance param- 
eter, Yp, is included in the likelihood to marginalize 


over an overall polarization efficiency. It scales CTE = 
p o ande eeu Xy. and we allow it to vary 


in a flat symmetric range centred at y, — 1.?! This pa- 
rameter is different from the frequency-dependent polar- 
ization efficiencies of C20, the frequency-dependent y,s 
have been combined in the CMB-only extraction step to 
retain the single overall polarization efficiency used here. 
The difference in efficiency amplitude across frequencies 
has also been accounted for in that step and therefore we 
here expect yp = 1. 

We provide both Fortran—90 and Python versions of 
the likelihood software, actpollite_dr4, with the 
DR4 data release on LAMBDA. 


6.2. Additional CMB data 


6.2.1. Large-scale polarization 


While our primary goal is to make a measurement 
of cosmological parameters that is independent of the 
Planck data, ACT by itself cannot constrain the op- 
tical depth to reionization, 7. ‘There is also evidence 
that the WMAP estimate of the large-scale polarization 
has dust contamination (Planck Collaboration XV 2014). 
We choose to include an estimate of 7 obtained in recent 
analyses combining Planck and WMAP large-scale po- 
larization measurements. We use as a baseline a Gaus- 
sian prior of r = 0.065 + 0.015, a conservative choice 
based on Planck-HFI, Planck-LFI (Planck Collaboration 
V 2018; Pagano et al. 2019) and Planck-LFI combined 
with WMAP estimates (Natale et al. 2020).?? 


21 We choose to work with symmetric ranges as this acts as an 
overall polarization calibration for the CMB bandpowers; we tested 
with the full multi-frequency likelihood that restricting polariza- 
tion efficiencies to be less than unity does not impact cosmological 
results. 

22 For reference, the optical depth estimated from Planck-HFI 
in the legacy release results is 7 = 0.054 + 0.007 (Planck Collab- 


6.2.2. WMAP data 


ACT has a minimum multipole of 600 in TT, and 350 in 
TE and EE, and so lacks data around the first two acous- 
tic peaks in TT and first full peak in TE/EE. Therefore 
in our nominal results (see §7) we include the full power 
spectrum information from WMAP using the public 9- 
year observations (Bennett et al. 2013) on large and in- 
termediate scales (2 < @ < 1200) in temperature and at 
intermediate scales in TE (24 < £ < 800) using the pub- 
lic WMAP 9-year likelihood software. We do not use the 
2 « l < 23 polarization likelihood?’ which is replaced by 
the T prior. 

ACT and WMAP only overlap on angular scales that 
are noise dominated in at least one of the two experi- 
ments and therefore we ignore correlations between the 
two datasets. We tested and confirmed this hypothesis 
with the same calculations described for Planck in the 
next subsection. The likelihood can then be well approx- 
imated as —2lnLacr — 2n LWMAP. 


6.2.3. Use of Planck data 


A goal of this paper is to give an estimate of the cosmo- 
logical parameters using measurements that are indepen- 
dent of Planck. However, we also study the consistency 
of the ACT and Planck data, and the impact of combin- 
ing the two datasets. For Planck we use the latest 2018 
TT/TE/EE high-/ plik-lite likelihood (Planck Col- 
laboration V 2018) (30 < £ < 2508 in TT; 30 < £ < 1996 
in TE/EE) in combination with the commander low- 
multipole likelihood in temperature (2 < £ < 29). In all 
cases we impose the same 7 prior used for ACT.?* 

To combine ACT with Planck we neglect the covari- 
ance between the two datasets and simply multiply the 
likelihoods. For this to be valid we truncate the multipole 


oration VI 2018), 7 — 0.0866 9 0505 in analyses using improved 
HFI maps (Pagano et al. 2019) and + = 0.051 + 0.006 from joint 
processing of LFI and HFI maps (Planck Collaboration 2020). The 
change in 7 due to the prior we use here causes a 1.66 shift in mean 
and 6096 widening of error in the correlated spectrum amplitude 
parameter and has very little impact on other ACDM parameters. 

23 Setting to false the use WMAP.lowl. pol flag. 

24 We note that because of this prior all Planck-derived results 
reported here will be slightly different from the Planck legacy re- 
sults in Planck Collaboration VI (2018). 
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range used for ACT, to keep the effect of double counting 
information below a chosen threshold. To choose these 
cuts we compute the Fisher-forecasted ACDM parame- 
ter errors derived by truncating ACT spectra below some 
minimum multipole range. We compute the ratio of er- 
rors when ignoring the AC'T-Planck covariance, to those 
where the full off-diagonal covariance is included.?? We 
then select the minimum multipole for which this ratio is 
below 1.05, i.e., where the impact of neglecting the ACT- 
Planck covariance is less than a 5% underestimation of 
errors. We find that multipoles below 1800 in TT should 
not be included for the ACT data in this case, but no 
cut in TE/EE is required. This @ range cut was found 
to also be valid for single-parameter ACDM extensions. 
To account for this we introduce a flag in the ACT likeli- 
hood to select the TT spectra at € > 1800 when combin- 
ing with Planck, and then set the joint likelihood to be 
—2InLacr — 2In£planck. In future analyses we expect to 
include the cross-correlation between ACT and Planck 
in a joint likelihood, to exploit the full angular range of 
both datasets. 


6.3. Parameter estimation methods 


To sample cosmological parameters we use the publicly 
available CosmoMC software (Lewis & Bridle 2002)?° and 
estimate the basic six ACDM cosmological parameters: 
the baryon and cold dark matter densities, Q,h? and 
Q..h?, an approximation to the angular scale of the acous- 
tic horizon at decoupling, Óyc, the reionization optical 
depth, 7, the amplitude and the scalar spectral index of 
primordial adiabatic density perturbations, A, and m;, 
both defined at a pivot scale ko — 0.05 Mpc^!. From 
these we derive and present constraints on the Hubble 
constant, Ho in km/s/Mpc, and on the amplitude of 
matter fluctuations on scales of 8h~'Mpc, og. We as- 
sume a single family of massive neutrinos with a total 
mass of 0.06eV and impose Big Bang Nucleosynthesis 
consistency relations to calculate the primordial Helium 
abundance. 

We also explore a set of extended models beyond 
ACDM that have a single parameter extension. A set 
of these are sensitive to the primary CMB temperature 
and polarization anisotropies: the effective number of rel- 
ativistic species, Neg, the primordial Helium abundance 
Yur, and the running of the spectral index dn,/dink. 
We also explore models that affect the degree of lensing 
in the power spectra: the total sum of neutrino masses 
X` m, and the spatial curvature Qg. 

MCMC chains are run with theory predictions com- 
puted using fmax = 6000 and stopped once the CosmoMC 
Gelman-Rubin convergence parameter, R — 1, reaches 
values smaller than 0.01. 


7. ACDM AND THE HUBBLE CONSTANT 
7.1. ACDM from ACT alone 


25 The off-diagonal covariance was estimated assuming that the 
CMB information contained in the ACT regions is a subset of that 
in Planck, ignoring differences due to data filtering or source mask- 
ing. 

26 The results were also reproduced with different codes for 
theory predictions and sampling, i.e., using CLASS+Cobaya, 
CAMB+Cobaya and CAMB+CosmoSIS (Lesgourgues 2011; Lewis 
et al. 2000; Zuntz et al. 2015; Torrado & Lewis 2020). 





Fic. 11.— Comparison of ACDM parameters (mean and 1o er- 
ror) estimated from different datasets: ACT alone (aslo shown in 
C20), ACT with the D220 prior, WMAP, Planck, ACT+WMAP 
and ACT+Planck (different colors). The light blue horizontal 
band highlights the lo measurement in our nominal case for 
ACT--WMAEP. Before including information from the first acous- 
tic peak, ACT prefers a best-fit value with higher ns and lower 
QpA? than WMAP and Planck, at 2.83-2.70 significance. We see 
excellent consistency between ACT+WMAP and Planck. 


We find that the ACT data are well fit by the ACDM 
model (as shown in C20), and our constraints on the six 
basic and two derived ACDM parameters are reported in 
Table 4 and shown in Fig. 11 and Fig. 12, compared to 
constraints from WMAP and from Planck. The polar- 
ization efficiency is tightly constrained around unity. 


21 


TABLE 4 
ACDM PARAMETERS WITH MARGINALIZED MEAN AND 68% CONFIDENCE LEVEL FROM ACT, ACT+WMAP, ACT+PLANCK AND PLANCK?®. 
THE BEST-FIT PARAMETERS FROM OUR BASELINE RESULT ARE SHOWN IN COLUMN 4. 





Parameter ACT ACT+WMAP ACT+WMAP ACT+Planck Planck 
best-fit 
Basic: 
1009h? 2.153 + 0.030 2.239 + 0.021 2.223 2.237 + 0.013 2.241 + 0.015 
1000. ^? 11.78 -E 0.38 12.00 + 0.26 12.12 11.97 + 0.13 11.97 + 0.14 
10f0Mc 104.225 + 0.071 104.170 + 0.067 104.180 104.110 + 0.029 104.094 + 0.031 
T 0.065 + 0.014 0.061 + 0.012 0.061 0.072 + 0.012 0.076 + 0.013 
Ns 1.008 + 0.015 0.9729 + 0.0061 0.9714 0.9691 + 0.0041 0.9668 + 0.0044 
In(1019 A,) 3.050 + 0.030 3.064 + 0.024 3.068 3.086 + 0.024 3.087 + 0.026 
Nuisance: 
Up 1.0008 + 0.0047 1.0033 + 0.0042 1.0028 1.0019 + 0.0046 — — 
Derived: 
Ho [km/s/Mpc] 67.9 1.5 67.6 1.1 67.1 67.53 + 0.56 67.51 + 0.61 
Og 0.824 + 0.016 0.822 + 0.012 0.828 0.8287 + 0.0099 0.8279 + 0.011 
Additional derived”: 
Qa 0.696 + 0.022 0.687 + 0.016 0.680 0.6871 + 0.0078 0.6867 + 0.0084 
to 13.832 + 0.047 13.772 + 0.039 13.786 13.791 + 0.021 13.791 + 0.025 
Sg 0.830 + 0.043 0.840 + 0.030 0.855 0.846 + 0.016 0.846 + 0.017 
1040* 104.252 + 0.071 104.189 + 0.067 104.199 104.128 + 0.029 104.112 + 0.031 
DÌ 13.972 + 0.091 13.860 + 0.058 13.839 13.879 + 0.027 13.878 + 0.028 
Tdra 148.6 + 1.0 147.06 + 0.63 146.89 147.18 + 0.29 147.14 + 0.30 





a Parameters estimated from Planck data alone (using T'T'TEEE and the same prior on 7 used for ACT) are reported for comparison; 
the official Planck results (e.g., Table 2 of Planck Collaboration VI 2018 for TT,TE,EE+lowE) have a stronger constraint on 7 and 


correlated parameters, including As and og. 


^ Dark Energy density, Q4, the age of the Universe in Gyr, to, a combination of og and the matter density, Sg = o0g,/Qm/0.3, the 
acoustic scale at the CMB last-scattering, 0*, the angular diameter distance to last-scattering in Gpc, D4, and the sound horizon at 


the end of the baryonic-drag epoch in Mpc, rarag. 


The residuals of the ACT data with respect to the 
best-fit ACDM model are shown in Fig. 13 for the TT, 
TE and EE spectra, together with the WMAP data in 
TT and TE at larger scales not measured by ACT. This 
best-fitting model has a y? for ACT of 279 for 254 (260 
bandpowers minus 6 freely-varying parameters) degrees 
of freedom corresponding to a PTE=0.13. The model 
provides a good fit to the TT, TE and EE spectra. 

The estimated ACDM parameters from ACT are 
broadly consistent with both the WMAP and Planck 
constraints, and the ACT data by themselves now put 
tighter constraints on all the ACDM parameters than 
the WMAP data do, significantly so for the peak angle, 
Oytc, from which parameters including the Hubble con- 
stant are derived. We find excellent consistency in the 
estimated Hubble constant between ACT and Planck, 
to be discussed further in $7.3. ACT measures many 
more acoustic peaks than WMAP and its current power 
spectrum signal-to-noise ratio is about 1/3 of Planck’s 
(see 8A). The ACT DR4 significant improvement in con- 
straining power compared to earlier ACT analyses (Siev- 
ers et al. 2013; Louis et al. 2017) comes from measuring 
the CMB peaks to higher precision in both temperature 
and polarization due to the increased sky area and higher 
sensitivity (see $A). As shown in C20, TT and TE now 
carry similar weight in constraining ACDM parameters. 

We find that the dominant difference between param- 
eters estimated from ACT compared to WMAP and 
Planck lies in the spectral index, ns, and the baryon 
density, Qh?. In the absence of data constraining the 
Sachs-Wolfe plateau and the first acoustic peak, there is 
a strong anti-correlation between these two parameters, 
as shown in Fig. 14. ‘This was also noted in Planck 
Collaboration V (2018) and Planck Collaboration VI 
(2018) when constraining cosmology from Planck data 
using only the smaller angular scales. A lower value of 


the baryon density, which damps the small-scale power 
spectrum, can be partially compensated by a higher value 
of ns, which tilts the spectrum to restore the small-scale 
power. Within that degeneracy direction, the ACT data 
prefer a model that has a lower value of the baryon den- 
sity, and higher spectral index, than WMAP or Planck. 

Figure 15 shows the best-fitting spectra for each 
dataset: the ACDM solution obtained from ACT by itself 
has a large enough tilt that it underestimates the ampli- 
tude of the first peak in TT compared to either WMAP 
or Planck. 'The amplitude of temperature fluctuations at 
( — 220 (D220, hereafter) inferred from the ACT best-fit 
model is 2.40 lower than the WMAP measurement. 

Since these large scales have been well measured, and 
consistently, by both WMAP and Planck, it is natu- 
ral to include this information as a prior. If we sim- 
ply impose a weak prior on the amplitude of the first 
peak, choosing a conservative Gaussian prior that dou- 
bles the WMAP uncertainty on the amplitude D220, 
with D220= 5750.6 + 71.1 uK?,?^ we find that this trun- 
cates the Q,h?—n, degeneracy and brings the parameters 
estimated from ACT in closer agreement with WMAP 
and Planck, as shown in the top panel of Fig. 14 with 
the ACT+D220 contours. 

A stronger prior is to include the full WMAP infor- 
mation. Before doing so we quantify the consistency of 
ACDM parameters estimated individually from ACT and 
WMAP data. As in e.g., Addison et al. (2016) we use 
a quadratic estimator to test whether the parameter dif- 
ferences are consistent with zero, computing 


\? = (mean; —meanz)* (Cov;+Cov2)~ *(mean;—meang), 


14) 


27 This is the value from the re-estimate of WMAP results ob- 
tained with the 7 prior used here. We note that D220 from Planck 
yields the same result. 
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Fic. 12.— Constraints on the six basic and two derived ACDM parameters from ACT, WMAP and Planck using a common 7 prior. 
The diagonal shows the 1-dimensional posterior distributions and the contours show the 68% and 95% confidence regions. Parameters 
from ACT are consistent with WMAP (Planck) to within 2.3 (2.7)o; the largest difference is in the Q,h?-ns plane. The Hubble constant 


agrees to within lo. 


where mean; and Cov; are the mean and covariance of 
parameters estimated from each experiment, 7. ‘To ac- 
count for the common information coming from the T 
prior we compute the statistic in a 5-dimensional param- 
eter space using QyA?, Q.h?, Ouc, ns and A, exp(—2r), 
and we determine significance assuming that the param- 
eter posteriors are Gaussian. We find that the ACT and 
WMAP ACDM parameters agree at the 2.30 significance 
level. Repeating the same test for ACT and Planck cos- 
mologies, we find a slightly larger difference, at the 2.40 


level, or at the 2.70 level if we remove common modes 
between the two datasets (i.e., by excluding ACT TT at 
£ < 1800). 

The WMAP and Planck best-fit models prefer a some- 
what higher (lower) amplitude of the spectrum in TE 
(EE) compared to ACT, in addition to the higher first 
peak in TT as shown in Fig. 15. This difference cannot 
be explained by calibration or by our current model for 
instrumental beam leakage. The bottom panel of Fig. 14 
shows that an artificial 5% change in the TE calibration, 
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Fic. 13.— Panel 1, 3, 5: ACT CMB TT, TE and EE power spectra for co-added deep and wide spectra complemented by WMAP large 
scale TT and TE measurements, together with the best-fitting ACDM theoretical model fit to either ACT or ACT--WMAP data. Panel 
2, 4, 6: residuals with respect to the best-fitting models, with x? — 279 (288) for ACT (ACT+WMAP). The addition of WMAP data 
pushes the best-fit model towards a lower (higher) value of ns (Qh?) than for ACT-alone, which produces overall negative TE residuals 
for ACT in the range 1200 < £ « 1700, raising the TE x? by 9. This feature may be a statistical fluctuation. 


relative to the TT calibration, could shift the ACT con- 
straints towards the WMAP and Planck measurements. 
However, at present we have no reason to introduce such 
a correction, and such a calibration of the polarization 
would be inconsistent with the EE data. We also note 
that a similar shift is not generated by any of the ACDM 
model extensions considered here. 

The mild tension between the datasets may simply be 
a statistical fluctuation, or could be an unaccounted-for 
systematic effect in the ACT data, or a hint of a preferred 
model beyond ACDM. We expect to be able to draw 


stronger conclusions from additional data from ACT that 
has already been collected. 


7.1.1. The smoothing of the acoustic peaks 


In the analysis of the Planck data, an unexpected fea- 
ture was a preference for an excess amount of gravita- 
tional lensing in the power spectrum compared to the 
amount expected for the given cosmological model, quan- 
tified by the parameter Az (Calabrese et al. 2008). Here 
Ay = 1 if the degree of lensing matches the model pre- 
diction; Az = 0 would correspond to no lensing. The 
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Fic. 14.— Top panel: The direction of the correlation between 
Q,h? and ns; changes when TT large-scale information is included 
(68% and 95% confidence levels shown): for WMAP and Planck 
a larger ns increases the 2nd to 1st peak height ratio, which can 
be compensated by increasing the baryon density. For ACT, and 
Planck limited to £ 1000, the increase in n4 boosts the small- 
scale power, now compensated by decreasing the baryon density. 
Including a prior on the first peak height (D220, grey contours with 
hatches) brings the ACT, WMAP and Planck best-fit models into 
close agreement. Bottom panel: The ACT constraints move along 
the Q5? — n5 degeneracy line if an artificial 596 calibration factor 
is applied just to the T'E data. 


effect of lensing is to smooth the features of the acoustic 
peaks, and add small-scale power. This lensing parame- 
ter was found by Planck to be higher than what would be 
expected at the 2.8c level (Ar, — 1.180 4- 0.065) (Planck 
Collaboration VI 2018), and higher than what was found 
using their reconstructed lensing signal. The significance 
of this excess was reduced with the addition of lensing 
potential data or, to a lesser extent, with different likeli- 
hood versions (Efstathiou & Gratton 2019). 

With the new ACT power spectra data, we find no 
deviation from the standard lensing effect predicted for 
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Fic. 15.— Comparison of the ACT, WMAP and Planck best- 
fitting cosmological model spectra and residuals, plotted against 
the ACT data. The models have consistent peak positions. The 
ACT-alone model under-estimates the first TT’ peak amplitude 
by 2.40, and has on average a lower (higher) TE (EE) best-fit 
spectrum than WMAP and Planck, at less than 30 level; also seen 
in Fig. 13. 


a ACDM model, with 


A, =1.01+0.11 ACT, (15) 


consistent with unity to within lo. The marginalized 
distribution for Ay is shown in Fig. 16, together with 
the constraints from the TT, TE and EE spectra alone. 
These each give consistent results, with the main con- 
straining power from both TT and TE. 

This result is consistent with a companion paper (Han 
et al. 2020) that studies the degree of lensing in the power 
spectrum for a subset of the ACT data, in the D56 region; 
by ‘delensing’ the spectrum to remove some of the impact 
of lensing directly from the maps, they find no deviation 
from the expected lensing signal. 


7.2. Baseline results: ACDM from ACT and WMAP 


With high confidence in the large scale measurements 
where WMAP and Planck are signal-dominated and 
agree well, we take the approach of supplementing the 
ACT data with the WMAP data for our baseline results. 
As an independent dataset to Planck we can assess how 
the preferred cosmological model compares. 
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Fic. 16.— Constraints on the peak smearing amplitude, Az, 
as measured by ACT, by the individual TT, TE and EE spectra 
from ACT, and ACT combined with WMAP. The ACT data are 
consistent with standard gravitational lensing: the parameter Ay 
artificially scales the lensing potential relative to the standard 
Ar, — 1 model prediction (shown as a vertical line). The Planck 
measurement is shown for comparison. 


We report constraints from ACT--WMAP in Table 4 
and Figs. 11 and 17. We find that ACT+WMAP to- 
gether provide stringent new constraints on the cos- 
mological parameters and are in excellent agreement 
with the Planck legacy results (Planck Collaboration VI 
2018). Each individual parameter differs by no more than 
1.1c, and similar degeneracy directions are seen between 
pairs of parameters. The scalar spectral index is less 
than unity at more than 4o significance, consistent with 
the Planck detection of a departure from scale invari- 
ance of the primordial fluctuations. As shown in the ap- 
pendix 8A, these two independently-measured datasets 
do not have the same balance of statistical power among 
different parts of the data: Planck has comparatively 
more weight from the TT spectrum, while ACT has most 
weight from the TE spectrum, and comparatively more 
constraining power from the EE spectrum. 

The best-fitting ACT--WMAP ACDM model has 
XAcrT7 288,7 with residuals shown in Fig. 13. The wors- 
ening of the ACT y? relative to ACT-alone (Ay? = 9: 
+2 in TT, +9 in TE and -2 in EE) comes, as described 
above, from the TT at large scales and the TE at inter- 
mediate multipoles (800 < £4 « 2000, see Fig. 13). As 
mentioned above, the significance of this effect is not yet 
large enough to draw conclusions about whether it is a 
statistical fluctuation, or a hint of a systematic effect, or 
potentially a preference for a different model. 


7.3. The Hubble constant 


By combining ACT and WMAP we obtain a CMB- 
derived estimate of the Hubble constant, within the 
ACDM model, that is derived from the same sky seen 


28 If we allocate to ACT the same degrees of freedom as for the 
AC'T-alone fit, i.e., assuming ACT and WMAP both contribute 
to constraining all the varying parameters, this corresponds to a 
PTE=0.07. If we instead assume that half of the ACDM parame- 
ters are constrained by WMAP the ACT PTE is 0.09. 
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by Planck but is independently measured.?? We find 
Ho = 67.6 + 1.1 km/s/Mpc ACT + WMAP. (16) 


This 1.6% inference of the local expansion rate of the 
universe is remarkably consistent with the Planck mea- 
surement: 


Ho = 67.5 + 0.6 km/s/Mpc Planck .?? (17) 


We also find a 2.2% precision measurement of Ho from 
ACT alone: 


Ho = 67.9+1.5 km/s/Mpc ACT, (18) 


or Ho — 68.2 + 1.5 km/s/Mpc with the D220 prior, and 
note that the ACT-derived estimate is stable to a number 
of analysis choices (see C20). Given that TE provides 
significant constraining power for Ho (see §A) we also 
note that our estimate is stable to tests run on TE; for 
example, applying an artificial 5% calibration factor to 
TE moves the A9 mean by only 0.15oc. 

Without the ACT data, the WMAP-derived esti- 
mate has an uncertainty twice as large (Hp = 70.0 + 
2.2 km/s/Mpc, Bennett et al. 2013; Hinshaw et al. 2013). 

This agreement between two independently-measured 
CMB datasets at a more comparable precision level, 
adds weight to the robustness of the CMB-inferred mea- 
surement. This new estimate of the Hubble constant, 
like Planck, is inconsistent with the z < 1 SHọoES 
Cepheids/supernovae based measurement (Riess et al. 
2019) at > 40 within the ACDM model (see Fig. 18). 
It also disagrees at 30 with the HOliCOW estimate 
from strong lenses (Wong et al. 2019) but agrees with 
the more recent TDCOSMO/HOLiCOW revised esti- 
mate including more flexible lens mass profile modeling 
and use of stellar kinematics data (Birrer et al. 2020). 
Like Planck, it is consistent with the TRGB-based mea- 
surement (Freedman et al. 2019), as well as with esti- 
mates derived from Baryon Acoustic Oscillations from 
galaxy clustering in combination with other non-CMB 
probes including Big Bang Nucleosynthesis and lensing 
(e.g., Aubourg et al. 2015; Addison et al. 2018; Abbott 
et al. 2018a; Cuceu et al. 2019). 

ACT's measurement of Hy comes from measuring the 
angular scale of the peaks in the temperature and po- 
larization power spectra, coupled with the peak heights 
which constrain the cold dark matter density and baryon 
density. An explanation for how this constraint works is 
given in Knox & Millea (2020). To summarize, the peak 
scale is set by the ratio of the sound horizon, rz, to the 
angular diameter distance to the last-scattering surface, 
D*,, with 6* = r3/D%,. The sound horizon is given by 
the integral of the sound speed, c,(t), over time before 
last scattering, t*, with 


LU dt. 
=f Sel, (19) 


It is set by the baryon and cold dark matter densities 
which affect the sound speed and the scale factor a(t) 


22 The only information retained from Planck is the conservative 
estimate of T. 

30 This is our estimate using the same 7 prior as ACT; Planck 
Collaboration VI 2018 finds 67.3+0.6 using their ‘lowE’ likelilhood 
to constrain T. 
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Fic. 17.— Parameter constraints as in Fig. 12 for ACT+WMAP versus Planck ACDM and derived parameters. There is close agreement 
between these two independently-measured CMB datasets, including the inferred Hubble constant. ACT-+Planck constraints are also 


shown for comparison. 


before recombination, and are measured via the relative 
peak heights in the power spectra. Figure 19 shows how 
the sound horizon is anti-correlated with the total matter 
density.?! 

By measuring the peak scale, 0*, and given the sound 
horizon, one then finds the angular diameter distance 


31 Following Knox & Millea 2020 we show Tdrag rather than r$, 
i.e., the sound horizon at the time of the baryon-drag epoch rather 
than the sound horizon at the time of CMB last-scattering. The 
two horizons differ by about 3 Mpc and the same trend is seen for 
both. 


to last scattering. In a flat universe, this is given as a 
function of redshift, z, by 


* 


E Ü ?* dz’ 
DiG)7 gv, E) s 


where E(z) = (Q,(1 + z)* -- G4(1 — z)? 4 Q4)9? is a 
function of radiation, matter and dark energy densities, 
and c is the speed of light. For given densities, this then 
allows an estimate of the Hubble constant: a universe 
with a smaller distance to the surface of last scatter will 
have a larger Hubble constant today. Since the densities 


© CMB © CMB (This work) Æ+ SNia 


SNla-TRGB m 2019) 


: ete cebtielg, (Riess 2019) 


‘ACT DR4+WMAP (Baseline) 
ACT DR4 
Planck PR3 (Planck VI 2018) 
ACT DR2+WMAP (Calabrese:2017) 


WMAP igen 2013) : 
| Pre-ACT DR4 





66 68 70 72 74 76 
Ho [km/s/Mpc] 


Fic. 18.— The Hubble constant estimated from ACT and 
ACT+WMAP (this work) is in excellent agreement with the mea- 
surement from Planck (Planck Collaboration VI 2018). The con- 
straints are compared to the WMAP-alone estimate (Hinshaw et al. 
2013) and combination of WMAP with previous ACT temperature 
data (Calabrese et al. 2017) to show the additional information 
coming from ACT DRA. Within the ACDM model these measure- 
ments agree with local measures using distances calibrated with 
TRGB stars (Freedman et al. 2019), but disagree with those cali- 
brated using Cepheid variable stars (Riess et al. 2019). 
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Fic. 19.— Correlation between the Hubble constant and the 
sound horizon at the baryon-drag epoch, Ho — rarag, showing 
the CMB-SHoES discrepancy (as in Knox & Millea 2020). Col- 
ored contours are the CMB 68 and 95% constraints from ACT, 
ACT+WMAP and Planck; the horizontal band is the Cepheids- 
based measurement from Riess et al. (2019). To reach the higher 
Ho measured by SHoES, the CMB needs a lower value of the matter 
density, Qm, which increases the sound horizon, as demonstrated 
for ACT-alone. 
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are not perfectly measured we see this interplay in the 
parameter correlations, highlighted in Fig. 19: a lower 
matter density would increase the distance to recombi- 
nation, so can be compensated with a higher Hubble con- 
stant to decrease the distance. An example is given in 
the Appendix B to show how a ACDM universe with 
Ho = 74 km/s/Mpc gives a poor fit to the data: the 
matter density is reduced to better fit the peak angles, 
but the peak heights cannot be sufficiently adjusted by 
varying the other parameters. 


8. EXTENSIONS TO ACDM 


The new ACT data provide us with the opportunity 
to explore a number of physically-motivated single pa- 
rameter extensions beyond the basic ACDM model. In 
doing so we look in particular for models that might re- 
duce the mild 2.30 offset between the ACT and WMAP 
data, and the slightly stronger 2.70 difference between 
the ACT and Planck best-fitting models. 

A summary of the constraints on all of the extended 
model parameters is reported in Table 5 and shown in 
Fig. 20. We do not find any significant deviations from 
ACDM, and highlight our main findings below. Figure 20 
also shows how the basic six parameters move during ex- 
ploration of extended models due to parameter degenera- 
cies. We note that the Hubble constant is never driven 
to very high values. 

Like for ACDM, our nominal ACT constraints are esti- 
mated using the CMB-only foreground-marginalized like- 
lihood, described in C20. Recall that when extract- 
ing this CMB-only power we make use of the earlier 
ACT MBAC data which includes 220 GHz observations 
and thus helps anchor the foreground model and re- 
duces the associated uncertainties at high multipoles. 
Slightly different results would be obtained using the 
multi-frequency likelihood which does not include that 
MBAC data (an example of this is shown for the effec- 
tive number of relativistic species, Neg, in Fig. 22). 


8.1. Lensing-dependent parameters 


The usual assumption in ACDM is that the geome- 
try of the universe is flat, with no curvature. With the 
ACT data we find this to be a good assumption, with no 
evidence for spatial curvature from CMB alone: 


Q, — —0.00153-516 ACT + WMAP (21) 


as shown in Fig. 21. This constraint on the curvature 
comes from the lensing information in the ACT power 
spectra: without lensing the curvature and the distance 
to the CMB, given by a combination of 0,h?, Q,h?, and 
Ho, are degenerate. The fact that ACT sees no evidence 
for an excess, or lack of, lensing compared to the ACDM 
model, is reflected in the lack of preference for non-zero 
curvature. The Planck power spectrum data alone pre- 
fer a non-flat model at the 2-80 significance (Planck Col- 
laboration VI 2018) (see also Handley 2019; Di Valentino 
et al. 2019); this is connected to the deviation from unity 
of the expected lensing signal (Az) and its significance 
is reduced with different sky fractions and Planck like- 
lihood choices (Efstathiou & Gratton 2020). Our new 
measurement from ACT, coupled with the reconstructed 
lensing signal from Planck, lends additional support to 
the explanation that the preferred non-zero curvature in 
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TABLE 5 
BEYOND ACDM PARAMETERS WITH 68% CONFIDENCE LEVEL OR 95% UPPER LIMITS FROM ACT, ACT+WMAP, AND ACT+PLANCK. 


Parameter ACT ACT+WMAP ACT+Planck Planck? 

Qi —0.003*5074 — —0.001*0 010 —0.0187* 9-015 —0.087 5-020 
Xm, |eV] < Sell <1.2 < 0.54 < 0.37 

Neg 2.42 + 0.41 2.46 + 0.26 2.74+0.17 2.97 + 0.19 
dns/dink  0.069+0.029 0.0128+0.0081 0.0023 + 0.0063 —0.0067 + 0.0067 
YHE 0.211 +0.031 0.220+0.018 0.232 + 0.011 0.243 + 0.013 


? Planck alone results ((TIT'TEEE with the same 7 prior) are reported for reference. 
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Fic. 20.— Constraints on the six ACDM and extended parameters, for single parameter beyond-ACDM models. The dashed horizontal 
lines are centered on our nominal ACT+WMAP ACDM results and on the standard model value for the extended parameters. Different 
colors span different models; circles, squares and diamonds compare ACT with ACT+WMAP and ACT+ Planck. 


the Planck power spectrum is a statistical fluctuation. 
A tighter constraint on the curvature comes from com- 
bining the CMB data with measurements of the baryon 
acoustic oscillations scale from galaxy surveys (as done in 
e.g., Planck Collaboration VI 2018); since this measure- 
ment is dominated by the BAO data we do not repeat 
the constraint here. 

ACT is also sensitive to the sum of the neutrino masses 
primarily through the degree of lensing in the power spec- 
trum. The higher the neutrino mass, the more the am- 
plitude of structure will be suppressed, and the smaller 
the lensing signal will be. ‘The limit on the neutrino 
mass sum from ACT and WMAP is X m, < 1.2eV 
at 95% confidence. This upper limit is higher than 
Planck from the power spectrum alone, and is also con- 
nected to ACT’s inferred lensing signal being slightly 
lower than Planck’s. When combining with other ex- 
ternal datasets (for example lensing potential measure- 
ments from Planck (Planck Collaboration VII 2018) and 
baryon acoustic oscillations from BOSS DR12 consen- 
sus, 6dFGS and SDSS MGS datasets (Alam et al. 2017; 
Beutler et al. 2011; Ross et al. 2015)) to break degenera- 
cles, in particular with the cold dark matter density, we 


find consistent constraints to Planck (see Fig. 21) with 
ACT+WMAP combined with external datasets giving 
Yo my <0.27eV at 95% confidence. 


8.2. Primordial parameters 


Here we investigate the running of the spectral index, 
the effective number of relativistic species, and the pri- 
mordial Helium abundance. The high-n, region allowed 
by the ACT data tilts the spectrum so that the ACT 
data by themselves show a mild preference for having 
less damping in the small-scale power spectrum than in 
the ACDM model. Since the dominant effect of all three 
of these extension parameters is to affect the degree of 
damping, this is reflected in the ACT data alone prefer- 
ring a number of species Neg that is less than 3.046, a 
running of the spectral index dn,/dlnk that is greater 
than zero at the ~2.40 level, or a primordial Helium 
abundance, Yur, that is less than the Big Bang Nucle- 
osynthesis prediction of 0.247 (Pisanti et al. 2008; Pitrou 
et al. 2018). 

The constraints are illustrated for Neg in Fig. 22: we 
find some impact from our choice of foreground model- 
ing (the best-fit model moves depending on whether the 
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FIG. 21.— Constraints on the CMB-derived curvature parameter (left) and the sum of neutrino masses (right). These limits are driven by 
the lensing information in the power spectrum, as well as the primary CMB signal for neutrino mass. With the new ACT data we find no 
evidence of deviation from flatness, supporting the interpretation that the Planck primary CMB constraint (shown in purple for reference) 
is a statistical fluctuation. The neutrino mass is poorly constrained from ACT alone without the large-scale information. ACT+WMAP 
combined with external datasets (Planck lensing and BAO) gives » m, « 0.27eV at 9596 confidence compared to < 0.19eV at 9596 
confidence from Planck with the same data combination and same 7 prior. The standard model value for a flat Universe and the expected 
neutrino mass in an inverted hierarchy are shown with vertical grey lines. 
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Fic. 22.— Constraints on the effective number of neutrino species (left) and the running of the spectral index (right). Increasing Neg 


enhances the damping and shifts the peak positions; the ACT data prefer Neg lower than the standard 3.046, with less damping, but is 
still consistent with the standard model. The ACT data also prefer a positive running of the index at 2.40, but including large scales from 
WMAP removes a preference for a scale-dependent index. The standard model values are shown with vertical grey lines. 


MBAC data is used to constrain the foreground param- 
eters); and there is a degeneracy with n, which will be 
limited if the large-scale data are included. For all three 
models, combining with WMAP or Planck results in a 
distribution that moves to a parameter space more con- 
sistent with ACDM predictions, at the 1.5-2.2c level for 
ACT+WMAP. We do not find evidence that any of these 
extensions significantly improve the goodness of fit of the 
combined ACT+WMAP dataset. 

The data strongly disfavor a model with excess neutri- 
nos, with Neg = 4 ruled out at 60, and Neg = 3.5 disfa- 
vored at 4c. Neg is correlated with the Hubble constant, 
since an enhanced Neg will increase the expansion rate, 
reducing the sound horizon. To hold fixed the CMB peak 
positions one can decrease the angular diameter distance 
to the CMB by increasing the Hubble constant. Allowing 


Neg to vary has the effect of lowering both Neg and the 
inferred Hubble constant along the Neg — Ho degener- 
acy line compared to the ACDM values, in the opposite 
direction from the possible tension. 


8.3. Alternative models 


Are any of the tensions in the cosmological data due to 
a missing component of the model or a new process? A 
series of forthcoming papers will investigate some more 
extended models. 

A process that decreases the sound horizon can resolve 
the tension with the Cepheid-derived Hubble constant. 
Two examples are introducing a self-interaction between 
neutrinos that delays free-streaming until shortly before 
recombination (Cyr-Racine & Sigurdson 2014; Kreisch 
et al. 2019), or introducing a period of early dark en- 
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ergy soon before recombination (Poulin et al. 2019; Hill 
et al. 2020). These models have distinct signatures in the 
small-scale CMB. 

We know little about the physics of the dark sector. 
CMB observations can constrain the cold dark matter in- 
teractions with baryons (Gluscevic & Boddy 2018). The 
signature is at small scales, beyond the resolution probed 
by Planck. In testing for deviations from fluctuations 
predicted by inflationary models, we can search not only 
for a running of the index, but also a model-independent 
deviation from a power law (e.g., Hlozek et al. 2012). 


9. SUMMARY AND OUTLOOK 


ACT has measured the temperature and polarization 
of almost half the sky at arcminute resolution. The data 
taken during 2013-2016 have been presented and anal- 
ysed here as ACT DR4. ACT has over four times the 
data comprising DR4 from three more seasons of obser- 
vations (2017-2019) which are in the process of analysis. 
During the 2020 season we are further expanding the fre- 
quency range taking more data with a lower-frequency 
channel. ACT is now producing maps over large areas 
of the sky, and we have developed optimal techniques for 
combining them with the Planck maps in order to pro- 
duce definitive sky maps with an extended angular range 
compared to Planck alone. 

By fitting the ACDM parameters to the ACT DR4 
deepest 15% of the sky we have demonstrated that the 
standard cosmological model remains consistent with ob- 
servations. Our measurement of the Hubble constant 
from ACT combined with WMAP agrees with that of 
Planck. We find no departures from Euclidean geometry 
and recover the expected amount of gravitational lens- 
ing in the power spectrum. The best-fit parameters from 
ACT-alone differ from those of Planck-alone by 2-380. 
This could be due to a statistical fluctuation, unknown 
systematic effects in our data or something new (noting 
that ACT and Planck measurements cover significantly 
different ranges of angular scales), and will be further 
explored with data already in-hand from ACT observa- 
tions between 2017 and 2019. As we analyze more data, 
our cosmological constraints from polarization will im- 
prove to the point where they become competitive with, 
and more powerful than, temperature, giving additional 
consistency checks on the standard cosmological model. 
Looking forward, ACT will continue to observe half the 
sky through 2021. 
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APPENDIX 
A. THE CONSTRAINING POWER OF ACT DR4 


In this Appendix we report signal-to-noise (S/N) estimates for the ACT, WMAP and Planck spectra, supplementing 
Fig. 20 in C20 that compares the ACT DR4 power spectrum error bars with those from Planck and WMAP. To 
assess the overall significance of the data we compute the x? of the TT, TE and EE spectra compared to null (i.e., 
how significant the measurement is when compared to no signal) and define S/N as its square root, given in Table 6. 
We estimate this for both Planck and ACT using their CMB-only likelihoods which account for the uncertainty due 
to foregrounds, calibration, polarization efficiency and beams. The equivalent for WMAP is estimated from the full 


WMAP likelihood. 


TABLE 6 
SIGNAL-TO-NOISE ESTIMATES FOR ACT, WMAP AND PLANCK CMB SPECTRA. 


S/N WMAP Planck ACT 
pr TE EE|TT Te. BE TEPPPIELDEBE|' PP Te. Ee. TP Te see 


273 23 12 | 606 182 236 646 175 76 112 185 







We find that: (i) the S/N for all measurements is greatest in TT; (ii) Planck has overall greater constraining power 
than ACT in TE and EE (even if limited to the same £Z range, i.e. / > 300); (iii) ACT has roughly comparable 
constraining power in EE and TE and comparable to WMAP in TT; (iv) ACT’s current overall constraining power is 
about 1/3 that of Planck’s. 

What is not captured in the above is the constraining power of each spectrum as a function of multipole. We explore 
this by taking two ACDM cosmological models with the most different inferred Hubble constants that fit the WMAP 
data to within 2c of the best-fitting model. We then take the difference between these two models, and in Fig. 23 
show the difference normalized by the error on the spectra measured by Planck and ACT. This then gives a relative 
measure of the S/N for estimating the Hubble constant, the parameter that we are particularly interested in, as a 


function of multipole. 


| Dou io ec ae posed cetus betae cocci me Ed 
|(—— ACT TT | | | | | |— ACT TE 


IAD,|/c(D;) 





500 1000 1500 2000 2500 3000 3500 4000 500 1000 | 1500 2000 2500 3000 3500 4000 
Multipole £ Multipole 7 


| —— Planck EE 
| |—— ACT EE 


|AD,|/a(D;) 





500 1000 1500 2000 2500 3000 3500 4000 
Multipole / 


Fic. 23.— ACT and Planck signal-to-noise contributing to the estimate of the Hubble constant as function of multipole. The relative 
constraining power of each experiment is shown with the difference between two theoretical models with widely-spaced Hubble constants 
(63 and 76 km/s/Mpc, both consistent with the WMAP data to within 2c) normalized by the power spectrum error bar. This highlights 
the new information in the ACT DR4 data and their complementary with Planck: the largest contribution from ACT comes from the TE 
spectrum in the range 1000 < £ < 2500 and the small scales in TT at £ > 2000 and EE at £ > 1000 provide additional information. 
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In comparing Planck to ACT, we find that the Planck data have the most constraining power at angular scales 
£ < 1800, with the ACT data constraining at smaller scales. In TE, the transition is at ~ 1000, and in EE at £ — 600. 
Figure 23 also shows the relative importance of TT, TE and EE. Despite the TT S/N itself being higher than TE and 
EEF, it is the TE that has higher S/N for ACT in terms of determining the Hubble constant, and other cosmological 
parameters. This is consistent with the ACDM parameter constraints reported in C20, and with previous studies in 
e.g., Galli et al. (2014). 


B. HOW THE ACT DATA DISFAVOR A HUBBLE CONSTANT Ho — 74 KM/S/MPC WITHIN ACDM. 


The latest Cepheid-based estimate of the Hubble constant gives Hp = 74.03 + 1.42 km/s/Mpc (Riess et al. 2019). 
Here we show how such a high value of Hp would compare with our data and our nominal ACDM fit. Figure 24 
compares the ACT+WMAP best-fitting ACDM model and its residuals compared to the ACT data, to the equivalent 
for a ACDM fit where Ho has been fixed to Hp = 74 km/s/Mpc. 

As described in $ 7.3, such a model will have a lower matter density to fit the well-measured peak scale. This high-Ho 
model also has a higher tilt and baryon density, attempting to compensate for effects of the changed matter density, 
resulting in a theory model that has more power in the first peaks and less power at small-scales than the best-fitting 
model. Overall this high- Ho fit degrades the ACT fit by Ay? = 55, with excess power in the residuals at TT and EE 
£ > 1000. 
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FIG. 24.— Comparison between the ACT--WMAP best-fit ACDM model and the ACT+WMAP best-fit model with Ho fixed to 
74 km/s/Mpc. The residuals show significant degradation of the fit to the ACT data. 


